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r  Historically,  surveillance  and  target  axqulaitlon  devices  have  operated 
In  the  visible  portion  of  the  spectrum  where  the  eye  was  the  dominant  sensor. 
The  eye  Is  still  s  critical  component  in  s  surveillance  or  target  acquisition 
system;  however,  the  eye  may  now  be  aided  by  an  infrared  (IR)  detection  and  im¬ 
age  diaplay  capability.  Also,  weapon#  systems  have  IR  seekers,  sensor  devices 
that  allow  the  weapon  to  detect  and  lock  on  to  a  target  that  has  a  thermal 
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20.  ABSTRACT  (Continued). 

^signature  significantly  different  from  its  surround.  The  thrust  in  the  United 
States  to  further  develop  IR  target  acquisition  and  surveillance  capabilities 
and  the  thrust  in  the  North  Atlantic  Treaty  Organization  (NATO)  community 
to  camouflage  or  conceal  critical  elements  at  fixed  Installations  from  such 
devices  have  focused  attention  on  the  need  to  better  understand  the  character 
of  the  thermal  IR  signatures  of  not  only  targets  but  also  their  surrounds.  To 
date,  considerably  more  effort  has  been  focused  on  targets  than  or  surrounds. 

The  Terrain  Surface  Temperature  Model  (TSTM)  presented  herein,  was 
developed  to  help  fill  the  void  in  the  understanding  of  thermal  IR  signatures 
of  natural  terrain  surfaces  and  of  some  cultural  features.  The  model  estimates 
temperatures  of  actual  or  hypothetical  material  systems  and  for  actual  or 
hypothetical  weather  conditions.  The  model  handles  sensible  heat  transfer, 
latent  heat  transfer,  the  impact  of  cloud  type  and  cover,  and  seasonal/ 
geothermal  heat  fluxes.  The  material  system  can  be  handled  as  a  multilayered 
medium  with  discrete  physical  and  thermal  properties  assigned  to  each  layer. 

This  report  documents  the  TSTM  by  presenting  a  discussion  of  the  mathe¬ 
matics  of  the  model,  a  discussion  of  the  computer  program  input  file  and  its 
operation,  and  the  results  of  a  sensitivity  analyses  end  limited  verification 
tests  conducted  with  the  model.  Typical  parameter  values  of  material  systems 
descriptors  used  in  the  model  are  included  in  an  appendix. 
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PREFACE 


The  study  reported  herein  was  conducted  by  personnel  of  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  from  1  Oct  1978  to  1  Jun 
1980.  The  study  was  done  under  Department  of  the  Army  Project  No. 
4A762730AT42,  Task  A4,  Terrain/Operations  Simulation,  Work  Unit  003, 
Electromagnetic  Target  Surround  Characteristics  in  Natural  Terrains,  and 
Department  of  the  Army  Project  No,  4A762719AT40,  Task  CO,  Theater  of  Op¬ 
erations  Construction,  Work  Unit  006,  Fixed  Installation  Camouflage 
Methods  and  Materials. 

The  study  was  conducted  under  the  general  supervision  of  Dr.  John 
Harrison,  Chief  of  the  Environmental  Laboratory  (EL),  and  Mr.  Bob  Benn, 
Chief  of  the  Environmental  Systems  Division  (ESD),  EL,  and  under  the 
direct  supervision  of  Dr.  Lewis  E.  Link,  Jr.,  Chief  of  the  Environmental 
Constraints  Group  (ECG),  EL,  The  development  of  the  mathematical  model 
presented  herein  was  accomplished  primarily  by  Dr.  Lee  Balick,  on  assign¬ 
ment  to  ECG  from  Colorado  State  University.  Assistance  was  received 
from  Messrs.  Randy  Scoggins  and  Curt  Gladen,  ECG,  and  Dr.  James  Solomon, 
Mississippi  State  University,  through  an  Army  Research  Office/Battelle 
Institute  grant  for  temporary  assistance.  This  report  was  prepared  pri¬ 
marily  by  Dr.  Balick  with  assistance  from  Dr.  Link,  Mr.  Scoggins,  and 
Dr.  Solomon. 

Directors  of  the  WES  during  the  conduct  of  the  study  were 
COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical  Direc¬ 
tor  was  Mr.  Fred  R.  Brown. 

This  report  should  be  cited  as  follows: 

Balick,  L.  K. ,  Link,  L.  E. ,  Scoggins,  R.  K. ,  and 
Solomon,  J.  L.  1981.  "Thermal  Modeling  of  Ter¬ 
rain  Surface  Elements,"  Technical  Report  BL-Bl-2, 
prepared  by  the  Environmental  Laboratory,  Water¬ 
ways  Experiment  Station,  in  collaboration  with 
Misslssipl  State  University,  for  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station,  CE, 

Vicksburg,  Miss. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


_ Multiply _ 

British  thermal  units 
(International  Table) 

degrees  (angle) 

feet 


_ By _  To  Obtain 

1055.056  joules 

0.01745329  radians 

0.3048 


metres 


THERMAL  MODELING  OF  TERRAIN  SURFACE  ELEMENTS 


PART  Is  INTRODUCTION 


Background 

1.  Camouflage  and  target  acquisition  have  opposing  functions,  one 
to  hide  and  the  other  to  seek.  They  have  a  common  denominator,  however, 
in  that  the  features  that  surround  the  target  to  be  camouflaged  or  iden¬ 
tified  (sometimes  called  the  background)  are  critical  in  both  the  hide 
and  seek  role.  An  equally  intimate  knowledge  is  needed  of  the  charac¬ 
teristics  of  both  the  target  and  the  background.  In  essence,  making 
something  match  the  background  and  discriminating  something  from  the 
background  are  inverse  problems  that  require  the  same  technology. 

2.  Historically,  surveillance  and  target  acquisition  devices  have 
operated  in  the  visual  wavelength  bands  where  the  eye  remains  a  dominant 
sensor.  In  the  past  decade,  thermal  Infrared  (IR)  technology  has  come 
of  age  providing  sensors  with  new  capabilities  for  target  acquisition 
and  presenting  a  new  threat  for  camouflage.  Optimizing  IK  sensors  for 
target  acquisition  or  optimizing  camouflage  measures  to  defeat  such 
se&j ors  requires  a  quantitative  understanding  of  the  thermal  IR  sig- 

’'ur.vs  of  both  targets  and  backgrounds.  Because  of  the  large  Impact  of 
such  factors  aa  air  temperature,  solar  Insolation,  wind,  and  cloud  con¬ 
ditions  thermal  signatures  can  vary  rapidly — an  additional  complication 
for  consistent  performance  of  either  target  acquisition  sensors  or 
camouflage  measures. 

3.  The  Army-Ulde  Ground  Target  Signature  Program  (AWGTSP)  is 
addressing  the  need  for  a  target-background  design  data  base  for  sensor 
design  snd  evaluation  through  a  three-part  program.  The  first  part 
deals  with  the  development  of  a  battlefield  IR  signature  model  that  will 
allow  extrapolations  of  target  and  background  signatures  to  varying 
environmental,  climatic,  and  seasonal  conditions  throughout  the  world. 
The  second  area  deala  with  updating  a  tactical  signature  library  to  fill 


critical  gaps  in  the  existing  empirical  signature  data  base.  The  third 
program  area  deals  with  susceptibility  analyses  and  is  designed  to  en¬ 
sure  that  vulnerability  of  all  Army  tactical  materiel  is  known  so  that 
effective  camouflage  can  be  brought  to  bear. 

4.  An  equally  important  problem  is  the  camouflage  of  key  elements 
at  fixed  installations,  a  responsibility  of  the  U.  S.  Army  Corps  of  Engi¬ 
neers,  for  which  similar  background  and  target  Information  is  needed. 

These  objectives  cannot  be  achieved  solely  on  the  basis  of  measured  signa¬ 
ture  data  or  performance  tests  in  the  field.  Modeling  and  simulation  of¬ 
fer  the  tissue  to  tie  measured  data  together  and  allow  extrapolations  to 
other  environmental  conditions,  sensor  types,  and  camouflage  materials. 

5.  Work  on  the  AWGTSP  has  resulted  in  considerable  progress  in 
computer  codes  for  predicting  the  performance  of  surveillance,  target 
acquisition,  and  terminal  homing  devices  and  target  signatures.  The 
target  models  have  ranged  from  simple  to  complex,  the  more  sophisticated 
approaches  using  combinatorial  geometry.  To  date,  targets  have  received 
considerably  more  attention  than  backgrounds.  Since  both  the  target  and 
target  surround  have  to  be  dealt  with  simultaneously  in  a  battlefield 
scenario,  a  compatible  and  equally  capable  target-surround  modeling  pro¬ 
cedure  is  needed. 

6.  A  previous  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  study  (Link  1979)  defined  a  research  approach  for  developing  a 
realistic  target-surround  signature  data  base  for  use  in  the  design  and 
evaluation  of  imaging  and  nonimaging  sensors  for  surveillance,  target 
acquisition,  and  terminal  homing  devices,  and  presented  preliminary  pro¬ 
cedures  for  predicting  terrain  surface  temperatures  of  typical  components 
of  battlefield  environmental  settings  with  special  emphasis  on  the  Fed¬ 
eral  Republic  of  Germany  (GE).  The  work  consisted  of  formulating  terrain 
and  climatic  data  bases  for  the  GE,  assembling  computer  models  for  pre¬ 
dicting  the  diurnal  temperatures  of  broad  classes  of  terrain  components 
as  a  function  of  climatic  conditions,  and  developing  a  matrix  for  inte¬ 
grating  the  GE  data  bases  and  the  temperature  models  to  provide  the 
initial  capability  to  predict  expected  temperature  ranges  for  terrain 
surface  features  under  a  variety  of  battlefield  scenarios. 
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7.  The  work  accomplished  in  the  above  referenced  study  provided 
an  initial  framework,  however  cursory,  for  estimating  thermal  regimes 
of  terrain  surface  features  for  specific  climatic  conditions.  Although 
the  data  base-model  framework  is  valid,  many  gaps,  both  data  and  ana¬ 
lytical,  were  evident.  These  gaps  and  the  perceived  capabilities  needed 
to  allow  complete  appraisal  of  the  performance  of  electro-optical  sys¬ 
tems  have  been  the  focus  of  additional  research. 

8.  The  work  that  followed  focused  primarily  on  an  enhancement  of 
the  mathematical  modeling  capabilities  illustrated  in  the  referenced 
report.  As  in  the  previous  study,  the  modeling  efforts  were  separated 
into  two  thrusts — models  for  vegetation  canopies  and  models  for  non- 
vegetative  surfaces  such  as  soil,  rock,  and  roadways.  The  vegetation 
canopy  modeling  has  been  accomplished  primarily  by  Colorado  State 
University  (CSU)  (Smith  et  al.  1980)  under  contract  to  the  WES.  The 
Thermal  Vegetation  Canopy  Model  developed  by  CSU  is  described  in  the 
reference  cited.  The  terrain  surface  (nonvegetative)  thermal  modeling 
work  was  done  primarily  inhouse  at  the  WES  and  is  the  subject  of  this 
report. 


Objective  and  Scope 

Objective 

9.  The  objective  of  the  work  presented  herein  was  to  generate  a 
capability  to  realistically  predict  the  temperature  (and  radiative  sig¬ 
nature)  histories  of  nonvegetative  natural  and  cultural  features  that 
commonly  comprise  the  backgrounds  to  targets.  The  ability  to  project 
background  temperatures/signatures  provides  a  means  to  examine 
temperature/signature  contrasts  that  occur  between  targets  end  back¬ 
ground  features  both  with  time  and  changing  weather  conditions.  This 

in  turn  provides  the  basic  Information  needed  to  examine  the  performance 
of  existing  or  proposed  target  acquisition  devices  and  the  effectiveness 
of  alternative  camouflage  measures. 

Terrain  Surface  Temperature  Model 

10.  The  Terrain  Surface  Temperature  Model  (TS1M)  was  developed 
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to  estimate  the  temperatures  of  actual  or  hypothatlcal  material  systems 
sod  for  actual  or  hypothetical  weather  conditions.  A  premium  was  placed 
on  simplicity  and  flexibility  with  respect  to  operational  constraints. 

In  short,  a  model  was  needed  that  considered  the  dominant  physical 
phenomena  that  influence  material  temperaturea  and  yet  was  reasonable  to 
use.  Some  of  the  characteristics  of  the  model  developed  to  meet  these 
requirements  are  as  follows: 

a.  Time  dependence  through  one  diurnal  cycle. 

b.  Air  temperature  aa  a  state  variable. 

c.  A  system  described  with  up  to  six  layers  of  uniform 
properties. 

d.  Precipitation  and  condensation  not  considered. 

e.  Spectral  characteristics  of  materials  not  considered. 

The  model  hendles  sensible  heat  transfer,  latent  heat  transfer,  the  im¬ 
pact  of  cloud  type  end  cover,  and  seasonal/ geothermal  heat  fluxes. 

Scone  of  report 

11.  The  TSTH  is  documented  in  this  report.  The  documentation  is 
comprised  of  a  discussion  of  the  mathematical  model  (Part  IX),  a  discus¬ 
sion  of  the  computer  program  input  deacrlptlon  and  lte  operation 

(Pert  XXX),  and  presentation  of  the  results  of  a  sensitivity  analysis 
and  limited  verification  testa  conducted  with  the  model  (Part  XV).  Xn 
Part  IX,  the  mathematical  framework  of  the  model  and  associated  assump¬ 
tions  are  presented,  followed  by  discussions  of  the  surface  and  bottom 
boundary  conditions  considered  and  the  mathematics  used  to  describe 
relevant  phenomena.  The  numerical  techniques  ut  d  to  solve  the  heat 
transfer  equation  ere  then  presented. 

12.  Xn  Part  XXX,  considerations  are  given  for  selecting  descrip¬ 
tors  of  material  systems  and  atmospheric  conditions  as  well  as  selecting 
options  for  computer  code  operation.  These  considerations  are  organised 
according  to  the  major  types  of  Inputs  to  tho  computer  code  since  vir¬ 
tually  all  material  system  and  atmospheric  condition  discriptors  and 
control  messages  for  program  options  are  included  in  the  input  data 
file.  The  sensitivity  analysis  in  Part  XV  is  subdivided  into  discus¬ 
sions  on  atmospheric  parameters,  material  system  parameters,  sad  control 
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parameters  for  the  numerical  solution  technique.  Model  verification 
tests  were  conducted  by  comparing  model  predicted  surface  temperature 
values  to  measurements  of  bare  soil  and  concrete  surface  temperature. 
Part  V  presents  conclusions  and  recommendations. 

13.  Twc  appendices  are  provided.  Appendix  A  presents  a  simpli¬ 
fied  flow  chart  for  the  computer  code,  and  Appendix  B  gives  typical 
values  for  many  of  the  material  system  descriptors  used  in  the  model. 


PART  JI:  MATHEMATICAL  MQDEL 
Introduction 


14.  A  general  schematic  of  the  TSTM  concept  is  shown  in  Figure  1. 
The  model  predicts  surface  temperatures  for  a  multilayered  (1-6  layers) 
system  by  determining  energy  transfer  in,  out,  and  through  the  system. 

A  basic  assumption  in  the  model  is  that  the  layers  and  the  environment 
above  them  are  horizontally  uniform;  i.e.,  the  most  significant  heat 
fluxes  are  vertical.  Thus  the  temperature  T  estimates  result  from 
solving  the  one-dimensional  heat  equation 


3T(z,t) 

3t 


3z 


subject  to  the  boundary  conditions 

n 

I  bit  -  0  at  x  -  0 
1-1 


£  Bu  -  0  at  z  -  b 
i-l 

where  the  observable  surface  is  z  -  b  ,  lower  surface  is  z  -  B  ,  o(z) 
is  the  diffusivity  and  both  bit  and  Bit  ,  i  »  1,2, ...,n  ,  denote 
heat  fluxes  at  time  t  . 

15.  Reliability  of  the  results  depends  upon  the  extent  to  which 
the  "hemal  characteristics  in  each  of  the  layers  can  be  approximated  by 
constant  values  and  also  is  strongly  dependent  upon  the  approximations 
of  bit  ,  i  -  l,2,...,n  ,  taking  place  at  the  surface  which  is  exposed 
to  environmental  heat  fluxes. 

Surface  Boundary  Conditions 


16.  The  boundary  condition  at  the  surface  is  estimated  with  a 
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heat  balance  equation  composed  of  six  energy  components:  (a)  insolation; 
(b)  radiant  energy  from  the  atmosphere  and  clouds;  (c)  radiant  energy 
emitted  by  the  top  surface;  (d)  sensible  heat  loss  (conduction  and  con¬ 
vection);  (e)  latent  heat  loss  (evaporation);  and  (f)  heat  conduction 
into  the  material.  A  value  for  each  of  these  components  is  calculated 
at  each  time  increment. 

Solar  energy  input 

17.  The  main  energy  input  to  the  surface  boundary  is  insolation. 
The  model  allows  the  option  for  measured  values  to  be  input  into  the 
program  or  (by  default)  the  values  can  be  estimated  by  a  series  of 
empirical  relations.  Adjustments  to  measured  data  to  account  for  differ¬ 
ences  of  surface  slope  and  cloud  conditions  must  be  made  external  to  the 
program. 

18.  If  the  option  of  having  the  model  estimate  insolation  values 
is  chosen,  a  series  of  calculations  is  executed.  All  solar  radiation  is 
considered  to  be  direct;  i.e.,  if  it  is  found  that  the  sun  is  not  shin¬ 
ing  directly  upon  the  surface,  then  it  is  assumed  that  there  is  no  solar 
radiation.  If  there  is  direct  radiation,  the  radiation  is  modified  to 
include  the  effects  of:  (a)  attenuation  due  to  atmospheric  gases  and  . 
water  vapor;  (b)  attenuation  due  to  cloud  cover;  (c)  orientation  of  the 
surface  to  the  sun;  and  (d)  reflection  at  the  surface. 

19.  In  order  to  determine  whether  direct  solar  radiation  occurs, 
the  sun's  position  relative  to  the  plane  of  the  surface  is  computed.  If 
there  is  no  insolation  at  the  surface,  the  calculation  stops  here.  The 
equations  for  sun-slope  geometry  for  a  horizontal  surface  are  those  of 
Small  (1977)  and  those  of  Sellers  (1965)  for  sloped  surfaces.  Chapter  3 
of  Sellers  (1965)  has  an  excellent  discussion  of  sun-slope  geometry 
calculations. 

20.  The  following  equation  is  taken  from  Khale  (1977)  and  is  used 
to  compute  insolation: 

S  •  (1  -  <0(1  -  A(u*,z)](0.349)Srt  cos  z 

»  5  O 

+  (1  -  ag)((l  “  «0)/(l  “  aoag) )(9. 651)So  coe  z 
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where 


S  -  the  solar  radiation  absorbed  at  the  ground  with  no 
cloud  cover 

a  -  surface  albedo 
g 

A(u*,z)  *  Mugge-Mbller  absorption  function,  equal  to  0.271(u* 

x 0.803 
sec  z) 

u*  ■  effective  water  vapor  content  of  atmosphere 

(0.349)S  ■  amount  of  solar  radiation  of  wavelength  greater  than 
0  0.9  pm 

Sq  ■  solar  radiation  incident  on  top  of  the  atmosphere 

z  -  zenith  angle  of  the  sun  as  a  function  of  time  of  the 
day  and  time  of  the  year 

a  *  atmospheric  albedo  for  Rayleigh  scattering,  equal  to 
°  0.085  -  0.247  log1Q  (<Ps/po)  cos  z] 

P8  ■  surface  pressure 

p  ■  1000  mb 
_o 

a  »  area  average  ground  albedo 

g 

(0.651)S  •  amount  of  solar  radiation  of  wavelength  less  than 
°  0.9  ym 

Excluding  clouds,  the  effective  water  vapor  content,  u*  ,  is  the  total 
precipitable  water  in  units  of  grams  per  square  centimetre.  Precipi- 
table  water  is  estimated  from  surface  air  temperature  and  relative 
humidity  with  the  following  equation  (Smith  1966): 

u*  ■  exp  (0.07074  4-  t] 

where  T.  is  the  dew  point  temperature  (degrees  Celsius)  and  r  ■ 
a 

-0.02290  from  April  to  June  and  t  -  0.02023  in  other  months. 

21.  Haurwltz's  (1948)  empirical  equation  is  employed  for  overcast 
sky  adjustments  with  respect  to  cloud  cover  and  type: 

CA  *  (a/94.4)  x  exp  l-tt  x  (b  -  0.059) J 

where  CA  la  the  cloud  adjustment  factor,  a  and  b  are  empirical 
coefficients  dependent  upon  cloud  type,  and  m  is  the  secant  of  the 
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solar  zenith  angle.  Values  for  a  and  b  are  available  for  the 
following  eight  cloud  genera:  cirrus,  cirrostratus,  altocumulus,  , 
altostratus,  stratocumulus,  stratus,  nimbostratus,  and  fog  (Table  1). 
Following  work  by  Pochop,  Shankiin,  and  Horner  (1968),  which  suggests 
interpolation  by  the  square  of  the  cloud  cover,  the  model  employs  the 
following  interpolatory  equation  for  intermediate  cloud  cover  situations: 

S  =  S  -  |[S  -  (S  x  CA)]  x  CC2 

c  a  \  a  a  1 

w>'ere  Sc  is  the  energy  reaching  the  upper  surface  and  CC  is  the  visual 
cloud  cover  in  tenths. 

22.  After  the  influence  of  clouds  is  estimated,  the  net  insola¬ 
tion  (S^)  is  adjusted  to  account  for  the  surface  slope  orientation  to 
compute  effeccive  in  .dent  net  insolation,  S  ,  as  follows: 

S  »  Sc  x  SF 

The  slope  factor  (SF)  i  defined  by: 

% 

SF  ■  cos  (z)  cos  (SI)  +  sin  (z)  sin  (Sl>  cos  (SAZ  -  S1AZ) 

where  z  is  the  solar  zenith  angle.  Si.  is  the  slope  of  the  surface 
(angle  between  slope  and  horizontal),  SAZ  is  the  solar  azimuth  angle 
and  S1AZ  is  the  azimuth  of  the  slope.  (Azimuths  are  measured  from  the 
south,  negatively  increasing  from  south  to  east  and  positively  increas¬ 
ing  from  south  to  west.) 

Atmospheric  and  cloud 
thermal  1R  energy  inputs 

23.  The  empirical  equation  used  to  estimate  atmospheric  IR  radi¬ 
ation  on  the  surface  II  is  the  Brunt  equation  (Sellers  1965): 

o 

U0  -  c«t‘  ^  +  b(e“-5)] 

where  c  is  the  emisblvlty  and  assumed  to  be  equal ,tr  1,  o  Is  the 
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Stephan-Boltzman  constant,  T  is  the  shelter  air  temperature  (degrees 
Kelvin),  e  is  the  water  vapor  pressure  (millibars)  and  b  and  c 

B 

are  empirical  constants.  The  values  of  Budyko  as  referenced  by  Sellers 
(1965)  are  used;  i.e.,  c  ■  0.61  and  b  -  0.050  .  The  value  of  e  is 

A 

obtained  from  Teten’s  equation  (Murry  1967) 

e  «RH  x  6.108  x  exp  (A  x  T  )/(T  +  273.15  -  B) 

B  B  B 

where  RH  is  the  decimal  relative  humidity,  A  »  17.269  and 
B  ■  35.86  .  The  cloud  contributions  to  thermal  IR  irradlance  I+fc 
are  treated  with  an  empirical  factor  adapted  from  Sellers  (1965): 

It.  -  1+  x  (l  +  CIR  x  CC2) 
t  o 

where  CIR  is  a  coefficient  dependent  upon  cloud  type.  Values  for  CIR 
can  be  found  in  Sellers  (1965)  or  Oke  (1978). 

Ground  radiative  emittance 

24.  The  surface  Is  treated  as  a  grey  body  emitter  such  that 

*♦  ’  V(T/ 

where  It  is  the  energy  radiated  from  the  surface,  c  is  the  cals- 

O 

aivlty  of  the  ground,  and  T  is  the  current  surface  temperature  as 

B 

predicted  by  the  model. 

Sensible  heat 

25.  The  conductive  and  convective  sensible  heat  transfer  H  is 
estimated  by  an  equation  following  Lamb  (1974)  or  Oke  (1978). 

h  -  or  *2«2  1£  &  scf 
H  -  -pCpte  i  H  SCF 

where 

(1.175(1  -  15Ri)0,75  Ri  <  0 

SCF  -  < (1  -  5R1)2  0  <  Ri  <  0.2 

(o  Ri  >  0.2 
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and  p  is  the  air  density,  is  the  specific  heat  of  dry  air  at 

constant  pressure,  k  is  von  Karman's  constant  (0.40),  z  is  the 
observation  height,  30/3z  and  3v/3z  are  the  partial  derivatives  of 
potential  temperature  and  wind  speed,  respectively,  with  respect  to 
height  z  .  Ri  denotes  the  Richardson  number.  Potential  temperature 
6  is  defined  by  the  relation 


e 


T 

a 


1000 \°* 286 
.  p  / 


where  T  and  p  are  the  air  temperature  and  pressure,  respectively. 
The  Richardson  number  is  defined  by: 

“  ■  (I  2)/®/ 

where  g  is  the  accelefatipn  due  to  gravity  and  6  is  the  average 
potential  temperature  between  the  surface  and  height  z  .  In  the 
model,  36/dz  and  3v/3z  are  approximated  by  first  order  differences, 
and  it  is  assumed  that  the  air  temperature  at  the  surface  equals  the 
surface  temperature  and  that  the  wind  velocity  is  zero  at  the  surface. 
Although  these  assumptions  at  the  surface  are  questionable,  they  elim¬ 
inate  the  need  for  and  difficulties  in  determining  aerodynamical  charac¬ 
teristics  of  the  surface.  This  formulation  was  selected  after  dis¬ 
cussions  with  Holbo,*  The  validity  of  the  assumptions  at  the  surface 
is  doubtful,  but  they  seem  preferable  to  the  selection  of  aerodynamic 
characteristics  (like  roughness)  and  attendant  assumptions  when  a  wide 
and  unpredictable  variety  of  surfaces  are  to  be  modeled. 

Latent  heat 

26.  The  following  equation  (Oke  1978,  Lamb  1974)  for  latent  heat 
exchange  E  is  used: 

E  -  -pLK2z2(W3q/8z)(3v/3z)SCF 


*  Personal  communication,  R.  Holbo,  Forest  Research  Laboratory,  Oregon 
State  University,  Corvallis,  Oreg. 
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-  V  H 


where  SCF  is  defined  above,  L  is  the  latent  heat  of  evaporation  (Hess 
1959),  q  is  the  specific  humidity,  and  W  is  the  decimal  relative 
saturation  of  the  top  surface.  In  calculating  a  value  for  3q/3z  ,  the 
air  at  the  ground  is  assumed  saturated.  The  value  of  q  at  z  *  Z 
represents  air  at  the  instrument  shelter;  i.e.,  3q/3z  represents  the 

potential  gradient  over  a  saturated  surface.  It  is  modified  by  the 
saturation  factor  W  to  account  for  nonsaturated  conditions.  As  with 
the  sensible  heat  loss  calculation,  E  is  set  equal  to  zero  if  the  air 
temperature  at  the  surface  is  less  than  or  equal  to  the  air  temperature 
at  the  shelter  height. 


Bottom  Boundary  Condition 

27.  The  bottom  boundary  condition  is  the  heat  flux  through  the 
bottom  of  the  lowest  layer  (top  layer  is  numbered  1,  second  layer  is  2, 
etc.  to  the  last  layer  (maximum  of  6>).  The  bottom  boundary  condition 
is  specified  by  one  of  the  following  three  options: 

a.  Option  1.  A  constant  temperature  at  the  bottom  boundary. 

h  ~ 

b.  Option  2.  A  constant  heat  flux  at  the  bottom  boundary. 

c.  Option  3.  A  constant  heat  flux  at  the  bottom  boundary 
and  an  additional  constant  temperature  radiating  surface 
below  the  bottom  boundary. 

28.  The  third  option  requires  additional  input:  (a)  bottom 
boundary  thermal  IR  emissivity;  (b)  bottom  boundary  geometric  shape 
factor;  (c)  under  surface  thermal  IR  emissivity;  (d)  under  surface 
geometric  shape  factor;  and  (e)  under  surface  temperature.  Geometric 
shape  factors  vary  from  0.0  to  1.0  and  are  related  to  the  emitting  and 
adsorbing  "efficiency"  of  the  surface.  High  efficiencies  result  from 
black  flat  surfaces  of  large  horizontal  dimensions.  Deviations  from 
this  configuration  suggest  lower  values  but,  unfortunately,  there  is 
little  basis  for  selecting  values  less  than  1.0.  Regardless  of  the 
option  chosen,  the  bottom  boundary  condition  is  kept  constant  in  time. 
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Numerical  Solution  „ 


29.  The  complicated  nonlinear  boundary  conditions  require  that 
the  heat  conduction  equation 


3T(z,t) 

3t 


3z 


be  solved  numerically.  In  this  equation,  u  is  the  diffusivity.  Each 
layer  is  assumed  to  be  homogeneous,  and  it  is  assumed  that  the  thermal 
characteristics  can  be  taken  to  be  constant;  specifically,  the  thermal 
conductivity  and  diffusivity  for  each  layer  are  assumed  to  be  a  constant. 
Solution  within  a  layer 

30.  Within  each  layer,  an  explicit  scheme  is  employed  to  solve 
the  one-dimensional  heat  equation.  In  particular,  given  the  temperature 
profile  at  time  t  ,  the  temperature  at  time  t  +  At  at  the  node  z  is 
given  by 


T(t  +  At,z)  •  T(t,z)  +  a(At/Az2) [T(t,z  +  Az)  -  2T(t,z)  +  T(t,z  -  Az)] 


where  At  is  the  time  increment  and  Az  denotes  the  spatial  Increment. 
It  should  be  noted  that  numerical  stability  requires  a  At/Az  <  1/2  . 
The  problem  of  numerical  stability  is  critical  for  thin  highly  conduc¬ 
tive  layers. 

Solution  at  the 
interface  of  two  layers 

31.  The  following  derivation  of  an  explicit  finite  difference 
scheme  to  handle  the  interface  between  layers  is  a  modification  of  that 
presented  in  Carnahan,  Luther,  and  Wilkes  (1964).  The  derivation  as¬ 
sumes  perfect  thermal  contact  at  the  interface,  l.e.,  continuity  of  the 
heat  flux  and  temperatures  at  the  Interfaces. 

32.  Let  layer  1  have  thermal  conductivity  k^  and  diffusivity 
ot^  and  layer  2  have  thermal  conductivity  and  diffusivity  kj  and 

02  ,  respectively. 
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LAYER  1  (kx. 


LAYER  2  (k2, 


a2) 


33.  Knowing  the  temperature  ,  T£  ,  and  1^  at  the  node 

points  i  -  1  ,  i  ,  and  i  +  1  ,  the  problem  is  to  calculate  the  new 
temperature  at  the  interface.  Employing  the  truncated  Taylor 
series,  is  approximated  by 


where  il  denotes  the  partial  derivative  in  layer  1  at  the  interface. 
Thus, 


Also,  the  first  order  approximation  to  3T/3t  is  given  by 


Since  3T/3t  -  c^O^/at2) 


one  obtains 
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'-Wl'  f* 


or 


34.  In  a  similar  fashion  for  layer  2,  one  obtains 


k 


2 


[TJ  -  TJ  + 


k  T 
11+1 

Az„ 


Vi 

Az„ 


Continuity  of  the  heat  flux  implies  that  Equation  1  equals  Equation 
thus,  after  simplification*  the  final  equation  used  to  calculate 


(2) 

2} 

is: 


Upper  boundary 

35.  The  new  or  updated  value  of  the  surface  temperature  T(t  + 
At,0)  is  calculated  by  solving  the  surface  heat  balance  equation  for 
the  surface  temperature  .  The  heat  balance  equation  is 

S  +  X*t-H-B-It  +  G»0  (3) 

where  G  denotes  tho  heat  flux  into  the  surface;  i.e.,  G  11  k(3T/3a) 
and  is  approximated  by  k(T^  -  T^/Az)  where  k  denotes  the  eonduc- 
tivlty  of  the  surface  layer  and  T^  denotes  the  temperature  at  the 
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present  time  for  the  first  node  point  below  the  surface. 
D  «  S  +  I+t  -  H  -  E  ,  the  heat  balance  equation  becomes 


•f  k 


0 


Letting 


or  upon  rewriting. 


kT, 


eaAz 


D_ 

eo 


0 


The  function  F  is  defined  by 


rdg) 


T  +  — -r — 
g  eaAz 


l  +  DAz 
coAz 


) 


(A) 


It  is  seen  that  the  updated  surface  temperature  is  a  root  of  F  .  The 
Newton-Raphson  algorithm  has  been  employed  to  locate  a  value  of  T 

8 

such  that  F  vanishes.  In  employing  the  Newton-Raphson  scheme,  the 

derivative  of  F  with  respect  to  T  is  needed: 

8 


dF(TR) 

dT 


8 


4T3  + 
8 


k 

eoAz 


dD/dT 
_ B. 

eo 


(5) 


Numerical  considerations  have  resulted  in  the  approximation  of  dD/dTg 
by  the  following  expression 


dD 

dT. 
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vv 

-AT 


(6) 


where  is  the  value  of  D  using  the  latest  estimate  of  T^  ,  Dq 
is  the  value  of  D  obtained  by  using  the  previous  estimate  of  T  ,  and 

O 

AT  denotes  the  change  In  temperature.  The  starting  value  for  the 
Newton-Raphson  scheme  is  taken  to  be  the  surface  temperature  at  the  pre¬ 
vious  time  step.  It  appears  that  three  to  five  iterations  yield  satis¬ 
factory  convergence  to  the  new  surface  temperature. 

Bottom  boundary 

36.  As  mentioned  in  paragraph  27,  the  bottom  boundary  condition 
is  the  heat  flux  through  the  bottom  of  the  lowest  layer  and  can  be 
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specified  with  one  of  three  options.  The  requirement  of  a  constant 
temperature  results  in  a  straightforward  boundary  condition. 

37.  For  options  2  and  3  (see  paragraph  27) ,  it  is  required 
that  the  following  equation  be  satisfied: 

R+-G-R+-D-0  (7) 


where  R+  denotes  the  radiative  energy  loss  through  the  bottom  boundary , 
G  denotes  the  heat  flux  into  the  lower  surface  and  is  given  by 
G  ■  k(3T/9x)  where  k  denotes  the  conductivity  of  the  bottom  layer, 

Rt  denotes  the  radiative  energy  from  the  constant  temperature  radiating 
surface  below  the  bottom  boundary,  and  D  is  the  constant  heat  flux  at 
the  bottom  boundary.  G  is  approximated  by  k(TR  -  T^)/6z  where 
is  the  temperature  of  the  bottom  surface  and  is  the  temperature  at 
the  first  node  point  above  the  bottom  surface. 

38.  The  following  equation  results  from  substituting  the  appro¬ 
priate  energy  components  into  Equation  7: 


Wb 


£R0bkRTR 


-  D 


(8) 


where  eB  denotes  the  bottom  boundary  thermal  IR  emissivity,  bRB 
denotes  the  bottom  geometric  shape  factor,  is  the  under  surface 
thermal  1R  emisaivity,  bRR  is  the  under  surface  geometric  shape  factor, 
and  T_  denotes  the  under  surface  temperature.  Equation  8  is  solved 

K 

by  employing  the  Newton-Raphson  iterative  scheme. 
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PART  III:  COMPUTER  PROGRAM  INPUT  DESCRIPTION  AND  OPERATION 


39.  The  discussion  presented  in  this  Part  is  intended  to  serve  as 
both  a  description  of  the  input  files,  options,  and  operation  of  the 
TSTM  computer  code  as  well  as  a  guide  for  users  of  the  model.  Since  all 
terrain  surface  feature  (material  system)  and  atmospheric  condition  de¬ 
scriptors  and  control  parameters  are  a  part  of  the  input  data  file,  the 
discussion  is  organized  around  the  major  types  of  inputs  to  the  model. 
Within  the  discussion  of  each  major  input  type,  the  input  file  format, 
parameters,  parameter  units,  and  relevant  control  options  are  presented. 
Guidance  is  blended  into  this  framework  for  data  input  selection  and 
consideration  of  the  available  control  options. 

AO.  The  input  data  are  read  by  the  program  as  a  single  data  file 
consisting  of  a  series  of  lines  of  data.  Each  line  of  data  starts  with 
a  line  number.  Following  the  line  number  are  the  necessary  data,  coded 
in  a  free-field  format.  The  program  prints  each  line  of  the  input  data 
file,  excluding  the  sequential  but  otherwise  arbitrary  line  number, 
together  with  variable  name,  units,  and  data  category  (1  through  8).  An 
example  is  presented  in  Figure  2a,  and  the  resulting  output  Is  given  in 
Figure  2b.  The  actual  program  coding  Is  not  presented  in  this  report. 
However,  Appendix  A  contains  a  program  flowchart  and  variable  definition 
list  for  use  by  readers  who  have  obtained  the  program  code.  Typical 
values  for  many  of  the  Input  constants  are  given  in  Appendix  B,  Govern¬ 
ment  agencies  can  acquire  a  copy  of  the  code  by  sending  a  request  to  the 
Commander  and  Director,  WES. 

Atmospheric  Condition  Inputs  and  Control  Options 
Atmospheric  constants 

Al.  In  addition  to  the  line  number,  the  first  line  of  input  data 
contains,  in  order,  atmospheric  pressure  (millibars),  wind  speed  (metres 
per  second),  cloud  type  index  (an  Integer  between  1  and  8),  and  meteoro¬ 
logical  Instrument  shelter  height  above  the  surface  (centimetres).  The 
shelter  height  Is  the  height  above  the  surface  at  which  air  temperature 
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and  wind  speed  are  measured.  Every  effort  mu^t  be  made  to  ensure  that 
the  air  temperature  and  wind  speed  data  reflect  the  conditions  at  the 
shelter  height  above  the  surface  (100-200  cm).  Means  of  adjusting  ob¬ 
servations  to  shelter  height  are  available;  see,  for  example.  Sellers 
(1965)  or  Hess  (1959).  Because  of  the  empirically  determined  coeffi¬ 
cients  in  the  relations  of  the  model,  values  for  shelter  height  should 
fall  between  100  and  200  cm.  Descriptions  of  cloud  genera  and  their 
corresponding  cloud  type  index  values  are  given  in  Table  1. 

Atmospheric  hourly  data 

42.  This  input  information  requires  24  lines  of  data.  Each  line, 
in  addition  to  the  line  number,  contains,  in  order,  the  24-hr  clock  time 
of  the  observation,  the  air  temperature  (degrees  Celsius),  relative 
humidity  (percent),  cloud  cover  (tenths,  0. 0-1,0),  wind  speed  (metres 
per  second),  and  total  insolation  (calories  per  centimetre  squared  per 
minute).  Recall  that  insolation  is  an  optional  input;  specifically,  the 
model  checks  the  input  data  to  determine  if  there  is  a  nonzero  value  for 
insolation  at  1200  hr.  If  a  zero  value  is  encountered,  the  program  then 
calculates  values  for  insolation  using  mathematical  relationships  men¬ 
tioned  in  Part  II  of  this  report.  It  is  clear  that  total  insolation  as 
an  accurate  measured  input  is  the  more  desirable  option  to  be  used  in 
the  model;  however,  effects  of  slope  changes  on  effective  insolation 
must  be  calculated  external  to  the  model  if  measured  values  of  insola¬ 
tion  are  input  into  the  program. 

Surface-Sun  Orientation  Specifications 

43.  In  calculating  the  seasonal  and  latitudinal  insolation 
changes  as  well  as  attenuation  due  to  atmospheric  gases  and  water  vapor, 
the  sun's  position  relative  to  the  surface  and  zenith  angle  is  needed. 

As  discussed  in  paragraphs  21  and  22,  the  orientation  of  the  sun  is  also 
employed  in  calculating  the  changes  resulting  from  cloud  cover  and  sur¬ 
face  slope. 

44.  For  this  line  of  data,  there  is  a  line  number  followed  by  the 
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angle  (degrees)*  of  surface  slope,  surface  azimuth  angle  (degrees), 

Julian  calendar  date  (integer  from  1  to  365) ,  and  latitude  (degrees) . 

45.  The  angle  of  surface  slope  is  merely  the  angle  between  the 
plane  of  the  surface  and  the  horizontal.  The  azimuth  is  the  direction 
of  the  vector  from  a  unit  zenith  vector  to  a  unit  normal  vector.  This 
model  uses  the  convention  that  south  is  0°,  increasing  positively  toward 
the  west  and  negatively  toward  the  east  (Sellers  1965).  The  Julian 
calendar  date  is  simply  the  number  of  days  that  have  lapsed  since 

31  December  of  the  previous  year. 

Heat  Flow  Calculation  Controls 

46.  This  line  of  input  data  contains  a  line  number,  followed  in 
order  by:  (a)  the  number  of  layers  (an  integer  from  one  to  six) ; 

(b)  the  number  of  24-hr  repetitions  to  be  run  (an  integer  greater  than 
or  equal  to  one) ;  (c)  the  time  increment  (minutes) ;  and  (d)  the  print 
frequency  (minutes). 

47.  The  number  of  24-hr  repetitions  gives  the  model  the  number  of 
daily  cycles  that  the  user  thinks  is  necessary  to  "stabilize"  the  initial 
temperature  profile.  If  the  initial  temperature  profile  is  considered 
accurate,  a  value  of  2  would  appear  reasonable;  whereas,  an  isothermal 
input  profile  may  require  four  cycles.  No  more  than  fiv'.  cycles  will 
ever  be  required  for  any  realistic  case.  Each  24-hr  cycle  is  run  with 
the  same  meteorological  conditions  and  parameters.  These  repetitions 
are  not  to  be  viewed  as  a  means  of  enhancing  the  numerical  calculations; 
i.e. ,  there  is  no  reason  to  think  that  by  increasing  the  number  of  24-hr 
cycles,  the  final  results  will  be  improved.  In  fact,  after  a  "stable" 
profile  has  been  obtained,  further  repetitions  could,  because  of  compu¬ 
tational  errors,  prove  detrimental. 

48.  The  time  increment  and  print  frequencies  are  entered  as 
integral  values  in  minutes.  In  selecting  a  time  increment,  the 

*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 
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stability  criterion  that  a  At/Az  <1/2  mast,  be  taken  into  considera¬ 
tion.  A  time  step  value  of  10  min  is  considered  to  be  a  useful  choice 
for  most  materials.  Thick  layers  of  good  insulating  materials  like  dry 
sandy  soil  may  be  well  represented  at  20-  or  30-min  steps,  but  thin 
layers  of  high  thermal  conductivity  may  require  time  steps  of  5  min  or 
less.  Tt  is  clear  that  the  print  frequency  twist  be  greater  than  or 
equal  to  the  time  increment  value. 

Material  System  Descriptors 
Initial  temperature  profile 

49.  This  set  of  data  lines  ia  composed  of  k  +  1  lines  where  k 
denotes  that  number  of  profile  points.  Specifically,  the  first  line 
contains  in  addition  to  the  line  number,  the  number  of  profile  points  to 
be  input.  The  next  k  lines  of  data  contain  the  depth  (centimetres) 
below  the  surface  and  the  temperature  of  the  material  at  that  depth. 

The  only  constraint  of  this  series  of  data  is  that  the  bottom  profile 
point  must  fall  within  the  bottom  layer.  As  discussed  in  paragraph  42, 
the  number  of  24-hr  repetitions  is  related  to  the  accuracy  of  the  tem¬ 
perature  profile  input  into  the  model. 

Top  surface  constants 

50.  This  line  of  data  contains  in  order,  following  the  ever  pres¬ 
ent  line  number:  (a)  the  thermal  IR  emisslvlty  (0.0  to  1.0)  of  the  top 
surface;  (b)  visible  wavelength  absorptivity  (0.0  to  1.0);  and  (c)  the 
moisture  content  (0.0  to  1,0). 

51.  Values  of  moisture  content  range  from  0.0  (dry)  to  1.0  (wet, 
saturated).  Because  of  the  complicated  nature  of  evaporation  for  non- 
saturated  surfaces  (Hess  1959),  the  user  should  be  cautioned  in  using 
Intermediate  values  for  this  variable.  In  fact,  the  user  might  consider 
using  something  such  as  interval  analysis  in  cases  where  intermediate 
values  are  used. 

Layer  specification 

52.  This  category  of  data  specifies  the  computation  parameters 
and  thermal  properties  fox  each  layer  of  the  material  to  be  modeled. 
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One  line  of  data  is  used  for  each  material  layer,  and  the  lines  are 
ordered  from  the  top  layer  to  the  bottom.  Each  line  of  data  has,  as 
usual,  a  line  number,  succeeded  in  order  by  the  following  information 
layer:  (a)  layer  thickness  (centimetres);  (b)  vertical  grid  increment 
(centrimetres) ;  (c)  thermal  diffusivity  (centimetres  squared  per  minute) 
and  (d)  thermal  conductivity  (calories  per  minute  per  centimetre  per 
degree  Kelvin) . 

53.  If  a  value  less  than  or  equal  to  zero  is  read  in  for  the 
vertical  grid  increment,  for  some  layer  the  program  will  use  a  default 
value  of  one-tenth  times  the  thickness  of  that  layer  for  the  vertical 
grid  increment.  Otherwise,  the  grid  increment  must  be  evenly  divisible 
into  the  layer  thickness. 

Bottom  boundary  specifications 

54.  As  indicated  in  paragraph  28  of  this  report,  the  user  has  the 
choice  of  one  of  three  options  for  the  bottom  boundary  conditions.  The 
initial  line  of  this  data  set  contains  a  value  for  the  Bottom  Boundary 
Index  (BBI),  which  informs  the  program  which  of  the  three  options  will 
be  used.  Specifically,  the  three  options  are  selected  and  defined  as 
follows : 

a.  Option  1.  BBI  <0  ;  e.g. ,  -1  indicates  that  a  constant 
heat  flux  through  the  bottom  layer  is  assumed.  If  this 
mode  is  chosen,  then  the  second  line  of  data  in  this  see 
contains  the  specified  heat  flux  value  (calories  per 
centimetre  squared  per  minute). 

1b.  Option  2.  BBI  -  0  ;  e.g.,  0  indicates  a  constant  temper¬ 
ature  at  the  bottom  boundary;  thus,  in  this  case,  the 
second  line  of  data  in  this  set  would  be  the  constant 
temperature  value  in  degrees  Celsius. 

c.  Option  3.  BBI  >  0  ;  e.g.,  1  indicates  the  option  of  a 
constant  heat  flux  at  the  bottom  boundary  and  an  addi¬ 
tional  constant  temperature  radiating  surface  below  the 
bottom  boundary. 

55.  Option  3  is  selected  if,  for  example,  there  is  air  space  be¬ 
low  the  bottom  boundary.  For  use  of  option  3,  the  second  data  line  re¬ 
quires  values  for  six  parameters:  (a)  the  constant  heat  flux  value 
(calories  per  centimetre  squared  per  minute);  (b)  the  bottom  boundary 
emlsalvlty  (0.0  to  1.0);  (c)  the  bottom  boundary  geometric  shape  factor 
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(usually  chosen  to  be  1) ;  (d)  the  lover  radiating  surface  emisslvity 
(0.0  to  1.0);  (e)  the  lover  radiating  surface  geometric  shape  factor; 
and  (f)  the  lover  surface  temperature  In  degrees  Celsius.  Emlssivities 
ape  the  familiar  grey  body  constants  that  range  from  0.0  to  1.0.  Geo¬ 
metric  shape  factors  are  intended  to  account  for  changes  of  thermal  IR 
emittance  and  absorption  efficiency  due  to  shape.  Large  flat  horizontal 
surfaces  have  high  values.  Small  area  sources  or  sinks  have  small 
values.  A  floor  in  a  large  room  vould  have  a  value  near  1.0,  but  a  hot 
motor  in  that  room  vould  have  a  smaller  value.  Unfortunately,  there  is 
little  guidance  for  an  a  priori  specification  of  values  for  the  geo¬ 
metric  shape  factors,  other  than  the  value  1.0. 
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PART  IV:  SENSITIVITY  ANALYSIS  AND  MODEL  VERIFICATION 


Sensitivity  Analysis 


Introduction 

56.  In  almost  any  application  of  the  TSTM,  uncertainties  in  one 
or  more  of  the  input  parameters  will  exist.  This  could  arise  from 
problems  with  instruments ,  uncertainty  of  the  materials  in  the  structure 
or  terrain,  or  a  number  of  other  possibilities.  Several  of  the  param¬ 
eters  are  difficult  to  obtain  accurately,  such  as  surface  saturation. 
Others,  such  as  cloud  cover,  require  decisions  by  the  user  which  can  be 
somewhat  subjective  in  nature.  All  this  will  usually  result  in  the  user 
dealing  with  a  range  of'  possible  temperature  predictions  rather  than  a 
simple  number. 

57.  The  range  of  uncertainty  in  model  output  can  be  established 
by  performing  a  sensitivity  analysis  on  the  variables  in  question.  The 
sensitivity  analysis  discussed  here  provides  a  starting  point  for  the 
user  to  detect  suspect  parameters  by  pinpointing  variables  to  which  the 
model  is  more  or  less  sensitive,  and  by  indicating  how  the  prediction 
might  change  for  variations  in  input  data.  This  analysis  was  performed 
using  a  standard  set  of  input  data,  while  varying  a  test  parameter  over 
a  range  of  values  characteristic  of  the  specific  variable.  Plots  were 
then  made  of  the  maximum  and  minimum  predicted  temperatures  versus  the 
value  of  the  test  parameter.  These  plots  are  contained  in  Figures  3 

to  31.  Those  curves  with  large  slopes  indicate  greater  model  sensitivity 
to  the  corresponding  parameter,  while  flatter  curves  indicate  little 
sensitivity.  A  summary  of  the  "standard  day"  used  here  is  given  in 
Figure  32.  Also  listed  are  the  material  parameters  for  the  three-layer 
system  modeled  in  this  analysis. 

58.  It  is  emphasized  that  this  sensitivity  analysis  was  performed 
for  only  one  set  of  materials  under  rather  artificial  conditions.  Also, 
since  the  analysis  of  all  possible  combinations  of  parameters  would  be 
too  enormous  to  be  practicable,  this  study  considered  univariate  changes 
only.  Therefore,  the  results  of  the  sensitivity  analysis  should  be  used 
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as  simply  a  guide  for  the  analysis  of  model  response  to  uncertainties 
in  the  input  data  for  each  caae  of  interest. 

Atmospheric  parameters 

59.  The  atmospheric  parameters  used  by  the  model  can  be  divided 
into  two  categories:  (a)  constant  and  (b)  time  dependent.  The  constants 
are  the  atmospheric  pressure,  cloud  type,  and  shelter  height.  These  con¬ 
stants  are  assumed  to  remain  unchanged  over  the  entire  24-hr  period  for 
which  the  model  is  run.  The  data  most  closely  corresponding  to  the  time 
of  interest  should  be  used.  The  model  exhibited  very  little  sensitivity 
to  pressure  and  shelter  height  for  all  values  of  saturation.  The  sen¬ 
sitivity  analysis  of  cloud  type  defined  three  groups  (see  Figure  11), 
types  1  and  2;  types  3,  4,  and  5;  and  types  6,  7,  and  8.  Within  each 
group  the  model  is  relatively  insensitive,  but  from  group  to  group  the 
changes  can  be  significant.  Of  course,  cloud  type  has  no  effect  in  the 
case  of  zero  cloud  cover. 

60.  The  time-dependent  atmospheric  parameters  are  the  air  tempera¬ 
ture,  relative  humidity,  cloud  cover,  and  wind  speed.  The  sensitivity 

of  these  variables  was  investigated  in  two  parts.  The  first  used  con¬ 
stants  for  all  hourly  values  in  each  run  of  the  model.  Secondly, 
diurnal  variations  were  introduced  into  the  input  time  series, 

61.  The  model  was  found  to  be  highly  sensitive  to  air  temperature, 
with  changes  in  the  maximum  and  minimum  predictions  that  were  of  the 
same  order  as  changes  in  the  air  temperature.  The  diurnal  variation  of 
the  air  temperature  produced  a  cycle  of  temperature  estimates  only 
slightly  damped  compared  to  the  amplitude  of  the  air  temperature  cycle. 
The  model *8  sensitivity  to  relative  humidity  occurred  primarily  in  the 
minimum  predicted  surface  temperatures.  The  maximum  values  change  only 
1.5*C  over  th<?  entire  range  of  relative  humidity,  while  the  minimum 
values  vary  5*C  (see  Figure  8).  Increasing  saturation  increases  maximum 
temperature  sensitivity  to  relative  humidity  variations. 

62.  For  cloud  cover,  predictions  were  made  for  all  cloud  types 
varying  cover  from  0  to  100  percent  for  each  type  except  8.  Type  8, 
fog,  has  only  0  or  100  percent  coverage.  The  model's  sensitivity  to 
cloud  cover  was  found  to  increase  with  an  increase  in  cloud  cover  and 
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in  the  index  values  of  the  cloud  type,  as  expected.  The  model  is  rela¬ 
tively  insensitive  to  any  cover  for  type  1,  but  temperatures  vary  more 
than  10°C  over  the  0  to  100  percent  range  for  type  7.  The  model  predic¬ 
tions  throughout  the  day  are  sensitive  to  cloud  cover  variations,  as 
seen  in  Figures  25  and  26. 

63.  The  sensitivity  of  the  model  to  wind  speed  is  primarily  in  the 
maximum  predicted  temperatures,  with  the  minimum  changing  less  than  1°C 
over  a  wide  range  of  wind  speeds.  However,  as  indicated  in  Figure  27 
and  28,  the  model  does  respond  to  change  of  wind  speed  with  time.  This 
is  also  true  for  cases  of  nonzero  saturation,  although  the  magnitude  of 
the  predicted  temperatures  decreases.  The  diurnal  variation  of  the  pre¬ 
dictions  indicates  that  the  model  is,  in  this  case,  insensitive  to  wind 
speed  variations  at  night  when  sensible  and  latent  heat  transfers  are 
minimal  or  zero. 

Site  characteristic  parameters 

64.  The  site  characteristic  parameters  consist  of  all  information 
supplied  to  the  program  describing  the  specific  structure  or  terrain  to 
be  modeled.  This  includes  data  on  the  geographic  location,  surface 
orientation,  time  of  year,  initial  conditions,  bottom  boundary  condi¬ 
tions,  and  material  properties.  The  model  was  not  tested  for  sensitiv¬ 
ity  to  surface  orientation,  geographic  location,  or  time  of  year  since 
these  are  exact  calculations  used  in  estimating  insolation.  However, 
users  who  must  estimate  insolation  in  cases  where  these  inputs  are  not 
exactly  known  should  consider  performing  a  brief  sensitivity  analysis  on 
these  parameters  for  their  particular  case.  Under  some  circumstances 
surface  temperatures  are  quite  sensitive  to  site  characteristics. 

65.  The  initial  temperature  profile  establishes  the  Initial  heat 
conditions  throughout  the  materials.  Relatively  small  changes  in  this 
profile  (thus,  material  heat  storage)  can  result  in  drastic  differences 
in  predicted  temperatures  (Figure  10).  Problems  with  this  parameter 
may  be  recognized  by  observing  the  predicted  temperatures  at  the  begin¬ 
ning  and  end  of  the  modeled  day.  Both  numbers  should  be  the  same,  if  the 
system  has  reached  equilibrium.  If  not,  the  user  may  increase  the  num¬ 
ber  of  24-hr  repetitions  made  before  printout  (see  Part  III).  The  best 
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solution  is,  of  course,  in  using  a  more  accurate  temperature  profile. 

66.  The  model  exhibited  little  sensitivity  to  the  bottom  boundary 
conditions  when  expressed  as  a  temperature,  varying  over  a  20°C  range. 
However,  the  effect  of  this  parameter  could  change  significantly  when 
imposed  at  smaller  depths  than  used  in  this  analysis  or  for  more  highly 
conductive  materials.  When  expressed  in  terms  of  a  net  heat  flux,  vari¬ 
ations  in  the  bottom  boundary  conditions  did  cause  small  changes  in  the 
predicted  temperature — increases  for  ingoing  flux  and  decreases  for  out¬ 
going.  Again,  this  effect  will  be  more  pronounced  for  shallower  bottom 
boundary  depths.  Wo  sensitivity  analysis  was  performed  for  the  case  of 
an  air  space  below  the  lower  boundary  (paragraph  28) ,  but  tests  were 
made  to  ensure  the  correctness  of  the  calculation. 

67.  The  material  parameters  consist  of  data  supplied  for  each 
layer  of  the  structure  of  terrain  and  information  needed  for  the  surface 
only.  Thickness,  thermal  diffusivity,  and  heat  conductivity  must  be 
specified  for  each  layer.  The  model’s  sensitivity  to  thickness  varia¬ 
tions  is  related  to  the  layer's  heat  conductivity,  with  higher  values 
resulting  in  less  change  with  thickness  than  lower  values.  For  constant 
layer  thickness  and  thermal  diffusivity,  increasing  heat  conductivity 
causes  a  narrowing  of  the  range  of  predicted  temperatures.  Variations 
in  thermal  diffusivity  (Figure  20)  produce  very  little  change  in  model 
predictions.  This  parameter  is  directly  involved  only  in  the  calcula¬ 
tion  of  the  numerical  stability  factor  described  in  Fart  II.  Values  for 
heat  conductivity  and  thermal  diffusivity  for  several  materials  and 
ground  cover  are  listed  in  Appendix  B. 

68.  The  surface  parameters  are  the  thermal  infrared  emissivity, 
solar  absorptivity,  and  'fractional  saturation  of  the  surface  material. 
The  model  exhibits  a  significant  sensitivity  to  emissivity  in  both  the 
minimum  and  maximum  predicted  temperatures,  with  both  decreasing  for 
higher  emisslvities.  Variations  in  solar  absorptivity  cause  large 
changes  in  the  predicted  maximum  temperatures,  which  increase  almost 
20*C  over  the  0  to  1  range  of  absorptivity.  The  minimum  temperature 
remained  relatively  constant,  however,  increasing  3°C  over  the  entire 
range  of  absorptivity  (see  Figure  23).  For  saturated  materials,  the 
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minimum  predicted  temperatures  change  very  little,  but  mayt-mim  values 
decrease  by  more  than  8°C  over  the  0  to  1  range  of  saturation.  Although 
not  investigated  in  this  analysis,  it  should  be  recognized  that  changes 
in  saturation  may  cause  variations  in  the  true  top  layer  heat  conduc¬ 
tivity,  thermal  diffusivity,  emissivity,  and  absorptivity.  No  provision 
for  automatically  modifying  these  parameters  for  effects  of  saturation 
is  included  in  the  program,  so  the  user  must  change  the  input  data  to 
account  for  these  effects. 

Numerical  procedures 
control  parameters 

69.  The  numerical  control  parameters  are  the  time  step,  the  grid 
spacing  (which  must  be  specified  independently  for  each  layer) ,  and  the 
number  of  24-hr  repetitions  to  be  performed  before  printout.  The  time 
step  and  grid  spacing  are  used  in  the  calculation  of  the  numerical 
stability  parameter  A  ,  as  shown  below: 


where 

■  thermal  diffusivity  for  i  layer 
At  *  time  step 

Az^  •  grid  spacing  for  i*"*1  layer 

In  order  for  numerical  stability  to  be  assured,  At  and  Az  must  be 

selected  so  that  A  is  less  than  or  equal  to  1/2.  For  a  large  diffu- 

2 

sivlty,  the  ratio  of  At  to  Az  must  be  small.  This  can  result  in 
very  small  time  steps  in  the  case  of  a  highly  conductive  thin  layer  be¬ 
cause  Az  mist  be  very  small.  Analysis  of  top  layer  thickness  indi¬ 
cates  that  thin,  highly  conductive  layers  may  usually  be  neglected  for 
multilayer  systems.  For  example,  in  the  case  of  a  2-mm- thick  metal 
plate,  a  temperature  gradient  of  only  0.5°C  is  created  across  the  layer 
(Bornemeier,  Bennett,  and  Horvath  1969).  Variation  of  At  and  Az  for 
constant  diffusivity  and  stability  factor  <1/2  shows  that  the  model  has 
essentially  no  sensitivity  to  these  parameters  provided  that  no  numer¬ 
ical  instability  is  encountered  because  At  or  Az  is  too  large. 
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Figure  24  shows  At  versus  Az  for  several  values  of  thermal 
diffusivity. 

70.  The  number  of  24-hr  repetitions  needed  is  directly  related  to 
the  accuracy  of  the  initial  temperature  profile.  The  analysis  of  the 
model's  sensitivity  to  this  parameter  indicates  that  when  the  initial 
temperature  profile  is  close  to  the  equilibrium  conditions,  predictions 
are  insensitive  to  repetitions  beyond  a  period  of  2  to  4  days.  Uncer¬ 
tainty  in  the  initial  profile  can  be  countered  to  some  extent  by  increas¬ 
ing  the  number  of  24-hr  repetitions. 

71.  In  Table  2  the  input  parameters  are  categorized  according  to 
the  degree  of  sensitivity  displayed  by  the  model.  The  model  exhibits 
little  sensitivity  to  some  parameters,  so  it  appears  that  relatively 
rough  approximations  can  be  used  for  these  with  confidence.  Above  all, 
any  user  should  recognize  that  this  model  is  not  expected  to  produce 
temperature  predictions  with  absolute  accuracy,  irrespective  of  the  cor¬ 
rectness  of  the  input  data.  There  are  simply  too  many  complex  processes 
in  nature  for  any  model  of  practical  usefulness  to  take  into  account. 

The  greatest  value  is  in  predicting  trends  for  such  uses  as  estimating 
contrast  between  two  dissimilar  objects  or  relative  changes  over  time. 
Take,  for  example,  the  case  of  contrast  prediction  between  a  structure 
and  the  surrounding  terrain.  Errors  from  poor  specification  of  atmo¬ 
spheric  parameters  will  be  common  to  both  features.  Here,  the  most 
important  parameters  would  be  those  which  are  different  for  the  structure 
and  terrain,  such  as  heat  conductivity  and  saturation.  The  user  must 
keep  considerations  such  as  these  in  mind  when  evaluating  the  uncer¬ 
tainties  of  model  estimates  for  a  specific  purpose. 

Model  Verification 


Technique 

72*  A  mathematical  model  is  only  useful  if  it  is  an  accurate 
abstraction  of  the  phenomena  modeled.  An  important  aspect  of  the  develop¬ 
ment  of  a  model,  therefore,  is  to  verify  that  the  model  provides  realis¬ 
tic  answers  for  the  conditions  or  bounds  within  which  it  is  to  be  applied. 
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One  method  of  model  verification  ia  to  compare,  values  predicted  by  the 
model  to  measured  values,  the  inputs  to  the  model  representing  the  con¬ 
ditions  under  which  the  measured  values  were  obtained.  This  is  a  rather 
simplistic  verification  because  only  the  final  model  output  is.  used  as  a 
yardstick  to  judge  the  ability  of  the  model  to  realistically  represent 
the  phenomena  abstracted  in  the  model.  This  method  also  dictates  that 
values  for  the  model  inputs  and  outputs  are  measured  very  carefully  and 
in  most  Instances  simultaneously.  As  already  pointed  out,  acquiring 
such  input  and  reference  output  values  is  often  a  very  difficult  task 
and  characteristically  such  data  are  seldom  readily  available  unless  a 
specific  effort  has  been  made  to  obtain  them.  In  many  instances  it  is 
difficult  to  accurately  measure  all  of  the  necessary  values  because  of 
a  lack  of  equipment  or  other  resources. 

73.  In  spite  of  the  shortcomings  mentioned  in  the  previous  para¬ 
graphs,  the  method  described  was  used  to  conduct  a  limited  verification 
of  the  TSTM.  Since  many  model  users  only  concern  themselves  with  model 
outputs,  usually  as  input  to  the  solution  of  some  broader  or  more  so¬ 
phisticated  problem,  it  does  provide  direct  and  visible  evidence  of 
model  applicability  and  accuracy. 

Test  data 

7A.  Two  terrain  surfaces  were  used  for  the  verification  test — a 
concrete  pad  and  a  bare  soil  patch  within  a  grass-covered  area,  both 
located  near  the  Environmental  Laboratory  Headquarters  building  at  WES. 
Weather  data  consisting  of  wind  speed,  wind  direction,  total  insolation, 
air  temperature,  relative  humidity,  and  precipitation  were  measured  at 
10-mln  intervals  for  the  period  28  July  80  to  5  August  80.  Surface  tem¬ 
perature  measurements  for  the  concrete  pad  and  bare  soil  surfaces  were 
also  made  at  10-min  intervals  coincident  with  the  weather  data.  The 
surface  soil  temperature  was  measured  with  a  thermistor  buried  approxi¬ 
mately  1  cm  below  the  soil  surface,  and  the  concrete  surface  temperature 
was  measured  with  a  thermistor  attached  to  the  surface  of  the  concrete 
pad. 

75.  All  weather  and  temperature  data  described  in  the  above  para¬ 
graph  were  recorded  on  Campbell  Scientific  Model  CR21  Microloggers. 
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All  weather  measurements  were  made  at  a  height  of  2  m  above  the  ground. 
Physical  entitles  such  as  concrete  pad  thickness,  soil  type,  soil  layer¬ 
ing,  initial  soil  temperature  profiles  (with  depth)  and  visual  reflec¬ 
tance  coefficients  of  the  concrete  and  soil  surfaces  were  measured. 
Values  for  heat  conductivity,  thermal  diffusivity,  and  absorptivity  were 
estimated  from  the  literature  (Link  1979).  Cloud  cover  and  cloud  type 
observations  were  made  periodically  during  daylight  hours  only  and  as¬ 
sumed  to  be  constant  for  both  days. 

76.  Weather  data  for  31  July  80  and  5  August  80  were  selected  for 
the  verification  test.  Two  days  were  selected  to  provide  some  measure 
of  repeatability.  Table  3  is  a  listing  of  the  weather  data.  Table  4  is 
a  listing  of  the  material  property  values  used  for  the  concrete  pad  and 
bare  soil.  The  concrete  pad  was  modeled  as  a  10-cm  layer  of  portland 
cement  concrete  over  a  15-cm  layer  of  sandy  soil.  The  bare  soil  patch 
was  modeled  as  a  single  50- cm  layer  of  uniform  silt. 

Model  predictions 

77.  The  TS1M  was  used  to  predict  surface  temperature  values  for 
the  bare  soil  and  concrete  pad  using  the  weather  data  for  5  August  80 
(Table  3)  and  for  the  bare  soil  using  the  weather  data  for  31  July  80. 
The  model-predicted  temperature  values  were  then  plotted  against  mea¬ 
sured  values  for  a  24-hr  cycle.  Figures  33,  34,  and  35  present  the 
predicted  and  measured  temperature  plots  for  the  concrete  pad,  bare 
soil  on  5  August  80,  and  bare  soil  on  31  July  80,  respectively.  Air 
temperature  is  Included  on  the  bare  soil  temperature  plots  (Figures  34 
and  35). 

Discussion 

78.  The  TSTM-predicted  temperatures  for  the  concrete  pad  (Fig¬ 
ure  33)  were  virtually  always  within  2*C  of  the  measured  values.  Per¬ 
haps  eveu  more  significant  was  the  close  match  of  the  predicted  and 
measured  values  for  the  peak  and  rising  and  falling  portions  of  the 
diurnal  cycle.  Throughout  the  nighttime  hours  the  predicted  values 
were  consistently  2*C  cooler  than  the  measured  values.  This  may  be 
due  to  such  things  as  cloud  cover  conditions  that  were  not  observed 
during  nighttime  hours,  or  an  insufficient  representation  of  sensible 
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heat  exchange  during  nighttime  hours.  Otherwise,  the  observed  tempera¬ 
ture  differentials  are  on  the  order  of  the  uncertainty  of  the  measured 
data. 

79.  The  predicted  and  measured  values  for  the  bare  soil  patch 
(Figures  34  and  35)  did  not  match  as  closely  as  those  for  the  concrete 
pad.  The  model  estimates  were  consistently  cooler  than  the  measured  val¬ 
ues  during  the  nighttime  hours  by  2  to  3°C.  Daytime  values  were  more 
closely  matched,  especially  for  the  31  July  80  diurnal  cycle.  In  both 
Figures  34  and  35  the  measured  values  lag  the  predicted  values;  this  is 
especially  evident  on  the  rising  portion  of  the  diurnal  cycle.  It  is 
hypothesized  that  this  lag  resulted  because  the  soil  temperature  measure¬ 
ments  were  made  at  about  1  cm  below  the  true  surface  (paragraph  74). 

There  is  characteristically  a  strong  temperature  gradient  with  depth 
below  the  surface,  and  the  temperature  fluctuations  below  the  surface 
lag  those  at  the  surface.  Thus,  the  actual  surface  temperatures  and  the 
model  predictions  of  those  temperatures  would  tend  to  respond  more 
rapidly  to  changes  in  the  weather  over  the  diurnal  cycle.  The  peak 
temperature  values  were  within  approximately  1  to  2*C  on  both  days, 
which  is  on  the  order  of  uncertainty  for  the  measured  values. 

80.  A  very  important  criteria  for  the  TSTM  la  the  capability  to 
predict  reasonable  temperature  contrast  values  for  different  terrain 
surface  scene  elements.  To  examine  the  ability  of  the  TSTM  to  generate 
contrast  values,  the  data  for  Figures  33  and  34  were  combined  and  plots 
of  the  temperature  difference  (contrast)  between  the  concrete  pad  and 
bare  soil  were  generated  for  both  the  predicted  and  measured  values. 
Figure  36  shows  a  comparison  of  the  predicted  and  measured  contrast 
values.  The  curves  in  Figure  36  very  closely  match,  indicating  that 
the  TSTM  did  a  good  job  of  predicting  contrast.  The  match  is  especially 
close  during  the  0  to  1600  hr  period.  The  relatively  larger  differences 
for  the  period  of  1600  to  2400  hr  can  be  explained  to  some  extent  by  the 
lag  of  soil  temperature  measurements  as  discussed  in  paragraph  79.  It 
is  important  that  the  model  was  able  to  represent  rapid  transient 
changes  of  contrast. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


81.  The  material  designer  or  evaluator  generally  chooses  to  work 
with  the  simplest  definition  of  an  operational  environment  that  seems 
practical.  In  the  case  of  surveillance,  target  acquisition  devices, 

or  camouflage  measures,  a  single  number  for  target-to-background  contrast 
is  often  the  criteria  used  for  design  or  evaluation.  Because  the  target 
and  the  background  are  often  both  dynamic,  the  contrast  value  is  also 
dynamic.  The  balance  between  simplicity  and  reality  must  be  carefully 
considered  and  every  effort  made  to  generate  the  capability  to  provide 
the  material  development  community  with  realistic  design  criteria  in  a 
cost-effective  manner.  The  AWGTSP  is  focused  on  this  broad  objective, 
and  the  mathematical  modeling  efforts  ongoing  such  as  that  reported 
herein  are  the  foundation  for  providing  the  tools  necessary  to  reach 
this  objective.  Realism  must  be  paramount  and  should  remain  the  major 
focus  of  data  base  development  and  modeling  efforts. 

82.  The  TSTM  is  a  step  forward  in  the  effort  to  realistically 
project  the  thermal  signatures  of  terrain  surface  (nonvegetatlve)  fea¬ 
tures  that  are  commonly  background  for  targets  or  items  to  be  camouflaged 
The  model  can  be  applied  with  either  actual  or  hypothetical  data  and  has 
been  shown  to  produce  realistic  outputs  over  a  range  of  atmospheric  and 
material  system  conditions.  Although  the  TSTM  is  a  one-dimensional  model 
it  can  be  used  to  examine  background  features  in  a  variety  of  situations. 
Difficulties  arise  in  attempting  to  model  vertical  walls  and  steeply 
sloped  surfaces.  Other  modeling  alternatives  are  being  examined  for 
these  situations.  Most  natural  nonvegetatlve  features  are  readily 
abstracted  to  layered  systems  with  horizontal  or  moderately  sloping 
surfaces. 

83.  The  major  advancements  offered  by  the  TSTM  ere  realistic  con¬ 
sideration  of  cloud  cover  and  cloud  type,  the  capability  to  simulate  the 
effects  of  fast  and/or  transient  changes  of  environmental  conditions, 
simplistic  inputs  describing  the  feature  and  atmospheric  conditions,  and 
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applicability  over  a  broad  spectrum  of  weather  and  seasonal  conditions 
and  material  characteristic  values.  The  major  difficulty  in  applying 
the  model  is  the  acquisition  of  accurate  values  for  material  descriptors 
such  as  heat  conductivity  and  thermal  diffusivity.  Some  materials  such 
as  portland  cement  concrete  can  vary  considerably  in  their  thermal  prop- 
erties  because  of  differences  in  the  characteristics  of  individual 
components. 

Recommendations 

84.  The  version  of  the  TSTM  presented  in  this  report  represents 
the  completion  of  a  phase  of  model  development.  Knowledge  gained  during 
this  phase  combined  with  an  increasing  awareness  of  its  short-  and  long¬ 
term  application  suggest  several  directions  for  continued  modeling  ad¬ 
vancement.  The  completion  or  a  model  development  phase  also  calls  for 
complementary  experimental  work  to  evaluate  the  present  model  and  to 
guide  future  effort. 

85.  The  model  contains  a  number  of  assumptions  and  simplifications 
which  potentially  limit  the  model’s  validity.  These  limits  need  to  be 
determined,  and  conditions  where  the  model  is  not  valid  must  be  defined. 
Theoretical  considerations  are  helpful,  but,  in  the  end,  the  model  must 
be  tested  against  carefully  obtained  experimental  data.  It  is  recom¬ 
mended  that  data  sets  be  obtained  and  analytical  effort  be  supplied  to 
provide  reasonable  demonstration  and  definition  of  model  validity.  (This 
is  no  trivial  task;  a  complete  experimental  validation  requires  data 
covering  all  combinations  of  all  controlling  variables  over  their  full 
range  of  values.  Accurate  measurements  are  not  easily  made.  Realistic 
time  and  effort  constraints  preclude  a  total  model  validation/ 
verification  effort.) 

86.  Layered  systems  are  common  on  the  earth's  surface,  making  the 
TSTM  a  widely  applicable  tool.  In  many  instances,  however,  the  terrain 
surface  is  partially  or  wholly  covered  with  a  grass  canopy.  While  vege¬ 
tation  canopy  models  are  available,  they  require  careful  specification 
of  intracanopy  meteorological  conditions,  canopy  structure,  and 


biophysical  characteristics.  The  general  can6py  signature  problem  is 
complex,  but  it  is  likely  that  a  simple  module  can  be  developed  for  the 
TSTM  that  can  handle  simple  herbaceous  canopies  under  moderate  environ¬ 
mental  conditions.  The  TSTM  could  then  be  applied  to  areas  such  as 
rangelands,  pastures  and  lawns.  It  is  recommended  that  a  simplistic 
herbaceous  canopy  module  for  the  TSTM  be  developed  and  tested. 

87.  Prediction  of  surface  temperatures  is  only  the  first  step  in 
the  prediction  of  thermal  IR  signatures.  Surface  temperatures  and 
spectral  emissivities  are  required  to  determine  spectral  radiance.  The 
effects  of  the  atmosphere  on  this  radiation  must  then  be  considered, 
followed  by  the  application  of  a  sensor  response  model.  Knowledge  of 
spectral  emissivity  of  terrain  materials  is  incomplete  and  unorganized. 

It  is  important  that  available  data  be  collected  and  that  the  remaining 
gaps  be  filled.  Atmospheric  transmission  modeling  work  is  currently  on¬ 
going  through  the  U.  S.  Army  Atmospheric  Sciences  Laboratory.  Their 
efforts  are  particularly  important  in  that  they  are  creating  an 
engineering- type  model  that  handles  the  complex  multiple  scattering 
phenomenon  that  is  significant  in  dusty  and  smoky  battlefield  environ¬ 
ments.  This  model  should  be  built  into  the  signature  prediction  process. 

88.  Consideration  should  be  given  to  the  development  of  an  Implicit 
scheme  for  the  solution  of  the  heat  transfer  equation.  Such  a  scheme 
would  be  unconditionally  stable,  thereby  expanding  model  applicability 
and  simplifying  model  operation. 

89.  The  TSTM  is  currently  one-dimensional.  While  this  configura¬ 
tion  can  handle  many  situations  (many  features)  that  occur  in  nature, 
there  is  a  class  of  features  that  may  require  a  two-  or  three-dimensional 
modeling  capability.  These  are  features  that  can  be  mistaken  for  targets, 
creating  false  alarms  or  false  target  acquisition.  Examples  might  in¬ 
clude  large  rocks,  cultural  features  with  dimensions  of  the  same  order 

of  magnitude  as  a  target,  or  small  Individual  facilities  that  are  to  be 
camouflaged.  It  is  recommended  that  research  be  conducted  to  investigate 
the  specific  need  for  a  multidimensional  modeling  capability  for  false 
target  features.  The  work  should  be  focused  on  the  magnitude  of  errora 
that  occur  from  applying  one-dimensional  models  for  these  features,  the 
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potential  Improvement  offered  by  multidimensional  models  development  of 
prototype  multidimensional  models,  and  more  verification  for  specific 
features. 
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Table  1 

Cloud  Genera  and  Cloud  Type  Indexes* 


Cloud  Genera 

Abbreviation 

Index 

Value 

Comments 

Cirrus 

Ci 

1 

High  clouds  composed  of  white  deli¬ 
cate  filaments,  patches  of  narrow 
bands,  elements  often  curved  or 
slanted  and  smaller  than  Cs, 
never  overcast  or  precipitating 

Clrrostratus 

Cs 

2 

High  clouds  appearing  as  whitish 
veil  usually  fibrous,  often  pro¬ 
duces  halo  phenomena,  thinner  than 
As,  does  not  appear  to  move, 
nonprecipitating 

Altocumulus 

Ac 

3 

Midlevel  clouds,  patches,  usually 
broken,  lee  wave  clouds,  elements 
smaller  than  Sc,  nonprecipitating 

Altostratu8 

As 

4 

Midlevel  grey  sheet  or  layer  of 
striated,  fibrous  or  uniform 
appearance,  large  horizontal 
extent;  thicker  than  Cs,  thinner 
than  Ns,  precipitation  generally 
light  and  continuous  (if  any) 

Stratocumulus 

Sc 

5 

Grey  and/or  whitish  layer  or  patch, 
nearly  always  has  dark  spots  and  is 
nonfibrous;  elements  larger  than 

Ac,  nonprecipitating 

Stratus 

St 

6 

Grey  rather  uniform  base,  patches 
ragged  if  present,  precipitation 
unusual  but  light  and  continuous 
if  present,  lower  and  more  uniform 
than  Sc,  less  dense  and  less  "wet" 
than  Ns 

NlmbostratU8 

Ns 

7 

Grey  often  dark,  diffuse,  large  hori¬ 
zontal  and  vertical  extent,  thicker 
than  As,  more  uniform  than  Sc,  of¬ 
ten  precipitating,  precipitation 
continuous 

Fog 

Fg 

8 

*  Cloud  genera;  Cumulus  (Cu),  Cirrocumulus  (Cc),  and  Cumulonimbus  (Cb) 
are  not  treated  here.  At  low  cloud  covers  (£0.3),  Cu  and  Cc  may  be 
approximated  with  Ac. 


ma 
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Table  2 

Relative  Model  Response  to  Variation  6f  Input  Parameters 


Very  Sensitive 


Moderately  Sensitive 


Very  Insensitive 


Air  temperature 


Relative  humidity 


Air  pressure 


Solar  absorption 

Thermal  emissivity 

Initial  temperature 
profile 

Saturation 

Cloud  cover  (types 
5,6,7) 

Top  layer  heat 
conductivity 


Shelter  height 

Wind  speed 

Cloud  cover  (3,4) 

Cloud  type:  (group  to 
group 
(1,2) 

(3,4,5) 

(6,7,8) 


Cloud  cover  (types 

1,2) 

Thermal  diffusivity 

Time  step* 

Grid  spacing* 

24-hr  repetitions 

Bottom  boundary 
flux** 


*  Not  sensitive,  provided  model  is  mimetic ally  stable. 

**  Not  sensitive  for  thick  systems  with  relatively  low  heat 
conductivity. 


Table  3 

Measured  Weather  Data  for  Model  Verification 


Tine 

Air 

Relative 

Cloud 

Wind 

Solar 

of 

Temperature 

Humidity 

Cover 

Speed 

Insolation 

day 

cC 

percent 

(0. 0-1.0) 

m/sec 

w/m2 

hr 

31  Jul/5  Aug  31  Jul/5  Aug 

31  Jul/5  Aug 

31  Jul/5  Aug 

31  Jul/5  Aug 

0 

23.5/24.8 

83.9/90.0 

0.3/0. 3 

0.8/0. 5 

0/0 

1 

22.3/24.5 

89.5/91.0 

0.3/0. 3 

0.5/0. 5 

0/0 

2 

21.7/24.1 

90.8/91.0 

0.3/0. 3 

0.5/0. 4 

0/0 

3 

21.6/23.9 

91.0/90.0 

0.3/0. 3 

0.4/0. 5 

0/0 

4 

21.5/23.6 

91.4/91.0 

0.3/0. 3 

0.5/0. 7 

0/0 

5 

21.6/23.2 

91.5/91.0 

0.3/0. 3 

0.6/0. 4 

0/0 

6 

21.0/23.2 

91.9/91.0 
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Table  4 

Material  Parameter  Values  for  Model  Verification 


Thickness,  cm 
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Thermal  diffusivity,  cm  /min 

Thermal  conductivity,  cal/min-cm-°K 

Longwave  emissivity 

Shortwave  absorptivity 

Surface  saturation 

Temperature  profile 
Depth,  cm;  temperature,  °C 

Bottom  boundary  temperature,  °C 
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CONCEPT  FOR  TERRAIN  SURFACE  TEMPERATURE  MODEL 


Figure  1.  Terrain  Surface  Temperature  Model 
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a.  Example  input  data  for  TSTM 
Figure  2.  Input  data  and  output  for  TSIW  (Continued) 
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b.  Output  for  TSTM 
Figure  2.  (Concluded) 
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BOTTOM  BOUNDARY  HEAT  FLUX,  CAL/CMJ-MIN 

Figure  7.  Changes  in  max-min  temperature  predictions  for  varia¬ 
tions  in  the  bottom  boundary  heat  flux 


Figure  8.  Changes  in  max-min  temperature  predictions  for 
variations  in  relative  humidity 
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INITIAL  TEMPERATURE  PROFILE, »C 


Figure  10.  Changes  in  max-min  temperature  predictions  for  varia¬ 
tions  in  the  initial  temperature  profile 


CLOUD  TYPE 

Figure  11.  Changt'3  in  max-min  temperature  predictions  for 
variations  in  cloud  type 
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CLOUD  COVER,  TYPE  5 

Figure  l6.  Changes  in  max-min  temperature  predictions  for  varia¬ 
tions  in  cloud  cover  using  cloud  type  5 


CLOUD  COVER.  TYPE  8 

Figure  17.  Changes  in  max-min  temperature  predictions  for  varia 
tions  in  cloud  cover  using  cloud  type  6 
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TOP  LAYER  THERMAL  DIFFUSIVITY,  CMJ/MIN 

Figure  20.  Changes  in  max-min  temperature  predictions  for 
variations  in  thermal  diffusivity 
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CASE  1  NO  CLOUDS  ALL  HOURS 

CASE  2  NO  CLOUDS  EXCEPT  OVERCAST  1100-1300  HR 

CASE  3  NO  CLOUDS  EXCEPT  OVERCAST  0800-1300  HR 
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Figure  25-  Model  predictions  for  different  diurnal  cloud 
cover  profiles*  base  cloud  cover  of  0  percent 
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CASE  4  WIND  SPEED  1  M/SEC  EXCEPT  10  M/SEC  0100  -  0600  HR 


Figure  i'?.  Model  predictions  for  different  diurnal 
wind  speed  profiles,  dry  surface 
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CASE  1  WIND  SPEED  1  M/SEC  ALL  HOURS 

CASE  2  WIND  SPEED  1  M/SEC  EXCEPT  10  M/SEC  1100  -  1300  HR 
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Figure  28.  Jtodol  predictions)  for  different  diurnal  wind 
speed  profiles),  saturated  surface 
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Figure  30.  Model  predictions  for  different  diurnal  relative 
humidity  p**ofilea,  dry  surface 
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Figure  31.  Model  predictions  for  different  diurnal  relative 
humidity  profiles,  saturated  surface 
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figure  32.  Standard  day  used  for  aensitivity  onolysio 
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Figure  33.  Predicted  and  measured  temperatures 
for  concrete  pad 

BARE  SOIL  (LOE  k>\  PATCH  IN  LAWN,  VICKSBURG.  MS,  5  AUG  80 
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Figure  31*.  Predicted  and  measured  temperatures 
for  bare  soil,  5  August  i960 
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Figure  35.  Predicted  and  measured  temperatures 
for  bare  aoil»  31  July  80 
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THERMAL  CONTRAST  (CONC,  TEMP.  -  SOIL  TEMP.) 
VICKSBURG  MS,  6  AUG  1980 
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Figure  36.  Predicted  and  measured  thermal  contrast 
for  concrete  pod  and  bare  soil 


APPENDIX  A:  SIMPLIFIED  FLOWCHART 
AND  VARIABLE  DEFINITIONS 


1.  The  logic  of  the  TSTM  computer  code -is  of  such  complexity 
(a  computer-generated  flowchart  requires  52  pages)  that  a  simplified 
flowchart  was  deemed  of  more  benefit  in  providing  the  user  an  understand¬ 
ing  of  the  organization  of  the  computer  code. 

2.  The  flowchart  presented  in  Figure  A1  avoids  the  many  "counters," 
such  as  those  found  in  the  location  of  node  points  and  layers,  time 
print  frequency,  and  interpolation  procedures,  in  the  program.  It  should 
be  noted  that  the  boundary  conditions  are  calculated  at  each  time  incre¬ 
ment;  thus,  each  of  the  energy  components  as  discussed  in  Part  II  is 
recalculated.  In  employing  the  Newton-Raphson  scheme  at  both  the  upper 
and  lower  boundaries,  the  energy  components  are  also  recalculated,  with 
the  last  estimate  of  the  appropriate  temperature,  for  each  iteration 
involved  in  obtaining  the  top  surface  temperature  and  the  temperature 

of  the  bottom  surface. 

3.  The  flowchart  is  followed  by  Table  Al,  which  gives  definitions 
of  the  major  variables  used  in  the  computer  code. 
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TO-CALCULATE-INSOLATION-ON-SLOPE  SURFACE 
“EMPIRICAL  RELATIONS  ARE  USED  TO  ESTIMATE 
INSOLATION  -  SEE  PART  II  OF  THIS  REPORT" 


TO  CALCULATE-TABLE-SLQPE-ANO-INTERCEPT 
“A  STRAIGHT  LINE  INTERPOLATION  IS  USEO 
BETWEEN  EACH  PROFILE  POINT  INPUT  INTO 
THE  MODEL" 


♦TO  CALCUlATE-SURFACE  ANO-LAYER  TEMPERATURE 


SET-UP-INITIAL-CONDITIONS 


“IN  EACH  LAYER,  TEMPERATURE  IS 
INTERPOLATED  FOR  EACH  NODE  IN 
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PRINT  OUTPUT-HEADING 
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Table  A1 

Variable  Definitions 


Variable 

ALPH(IX) 

APRM 

ATERM 

B 

BBB(J.I) 

BEP 

BX 

BPRM 

BTERM 

CLOUD 

DAY 

DECL 

DELT 

DIET 

DPRK 

DPRHO 

DPRH1 

DTERM 

ELF 

EPSN 

FACTA 

FACTO 

FACTE 


_ Definition _ _ 

THERMAL  DIFFUSIVITY  OF  LAYER  IX  IN  CM**2/MIN 

FACTE*TEMP**3  IN  CAL/MIN**2-CM**2 

ENERGY  CONTRIBUTED  BY  ATMOSPHERIC  IR  EMISSION 
CAL/CM**2-MIN 

HEAT  CONDUCTIVITY  OF  SURFACE  CAL/CM**2-MIN-K 

Y  INTERCEPT  OF  LINEAR  EQUATION,  USED  FOR  TABLE 
INTERPOLATION 

BOTTOM  SURFACE  EMISSIVITY 

BOTTOM  SURFACE  GEOMETRIC  SHAPE  IN  FRACTION  (0. 0-1.0) 

HEAT  CONDUCTIVITY  OF  BOTTOM  BOUNDARY  LAYER 

ENERGY  CONTRIBUTED  BY  INSOLATION  AFTER  ADJUSTMENT  USING 
SURFACE  ABSORPTIVITY.  IN  CAL/CM**2-MIN 

CLOUD  COVER  IN  FRACTION  OF  0. 1-1,0 

JULIAN  DAY  USED  IN  SOLVING  INSOLATION 

SOLAR  DECLINATION  ANGLE 

TIME  STEP  IN  HOURS 

DEPTH  IN  CM  OF  INITIAL  SOIL  PROFILE  AT  WHICH  CORRESPONDING 
SOIL  TEMPERATURE  IN  KELVIN  IS  INTERPOLATED 

HEAT  FLUX  IN  CAL/CM**2-MIN  BOTTOM  BOUNDARY  OR  TEMPERATURE 
IN  KELVIN  AT  BOTTOM  BOUNDARY 

TEMPERATURE  OF  BOTTOM  MATERIAL  IN  DEGREE  CELSIUS.  USED 
WHEN  LFLUXY  -  0 

HEAT  FLUX  THROUGH  BOTTOM  MATERIAL,  IN  CAL/CM**2-MIN, 

USED  WHEN  LFLUXY  NOT  EQUAL  0 

ENERGY  LOSS  DUE  TO  EVAPORATION 

LATITUDE  IN  RADIANS 

BMISSIVm  OF  SURFACE  MATERIAL 

SIGMA* EPSN 

FACTD-S IGMA*BK*8EP*TR**4  USED  IN  BOTTOM  BOUNDARY  CALCULA¬ 
TION  WHEN  THERE  IS  AIR  SPACE  BENEATH  THE  BOTTOM 

PACTE*S ICMA*RK*REF 

(Continued) 
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Table  A1  (Continued) 


Variable 

FACTH 

FK(IX) 

FMM(J,I) 

HEADER 

HTERM 

IEFSWT 

IEOF 

IMATL 

INTR(IX) 

IPFJNT 

IT1ME 

IX 

IY 

JMAX 

LAT 

LFLUXY 


L’ 

M 

MAX(J) 

NCLOUD 

NOMATL 

NPRNT 


_ Definition _ 

USED  IN  SOLVING  CONVECTION  TERM  (HTERM)  (1000.0/PRESS)** 
0.286 

HEAT  CONDUCTIVITY  OF  LAYER  IX  IN  CAL/MIN-CM-K 

SLOPE  OF  LINEAR  EQUATION,  USED  FOR  TABLE  INTERPOLATION 
OF  HOURLY  INPUT  DATA 

72  CHARACTER  INPUT  VARIABLE  USED  TO  PRINT  COMMENTS  ON 
OUTPUT 

ENERGY  LOSS  OR  GAIN  DUE  TO  CONVECTION  CAL/CM**2-MIN 

SWITCH  WHEN  «  0  WILL  PRINT  OUTPUT  ONLY  AT  SPECIFIED  TIME. 
IF  NOT  *  0  WILL  PRINT  OUTPUT  AT  EVERY  ITERATION 

SET  FROM  0  TO  1  WHEN  IEOF  IS  ENCOUNTERED,  USED  TO 
TERMINATE  PROGRAM. 

BACKWARD  COUNTER  OF  LAYERS.  STARTING  WITH  THE  NUMBER 
OF  LAYERS 

BEGINNING  SUB-LAYER  DEPTH  NUMBER  FOR  LAYER  NUMBER  IX 

BACKWARD  COUNTER  SET»NPRNT.  WHEN  EQUAL  TO  1  OUTPUT  IS 
PRINTED 

BACKWARD  COUNTER  INITIALIZED  AS  TOTAL  TIME  STEPS  IN  HOURS 

LAYER  NUMBER  STARTING  WITH  TOP  LAYER 

SUB-LAYER  DEPTH  NUMBER 

THE  TOTAL  NUMBER  OF  SUB-LAYERS 

LATITUDE  USED  IN  SOLVING  INSOLATION 

INPUT  BOTTOM  BOUNDARY  DATA  CONTROL  SWITCH.  IF  *  0,  THERE 
IS  NO  HEAT  FLUX  THROUGH  BOTTOM  OF  MATERIAL,  IF  NEGATIVE 
THERE  IS  NO  AIR  SPACE  BENEATH  BOTTOM  MATERIAL,  IF  POSITIVE 
THERE  IS  AIR  SPACE  BENEATH  BOTTOM  MATERIAL 

DUMMY  VARIABLE  TO  READ  LINE  NUMBER  FROM  INPUT  FILE 

SECANT  OF  SOLAR  ZENITH  ANGLE  IN  RADIANS 

NUMBER  OF  INPUT  TABLE  VALUES  USED  IN  TABLE  INTERPOLATION 
MODULE 

CLOUD  TYPE  INDEX  NUMBER  (1-8)  USED  IN  SOLVING  INSOLATION, 
INFRARED  EMISSION 

NUMBER  OF  MATERIAL  LAYERS  USED  IN  SOLVING  HEAT  FLOW 

NUMBER  OF  TIMES  OUTPUT  TIME  PRINT  FREQUENCY  IS  DIVISIBLE 
BY  TIME  STEPS.  USED  TO  DETERMINE  WHEN  TO  PRINT  OUTPUT 

(Continued) 
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Table  Al  (Continued) 


Variable 


Definition 


NTABL 
NX (IX) 

PRESS 

PTYHE 

REP 

RHOC(IX) 

RI 

RK 

RR(IX) 

SAZ 

SFRQ(IX) 

SICF 

SLOPE 

SMALLA 

SOLCHX 

SPEED 

ST0R(1,IY) 

STOR(2,IY) 

ST0R(3,IY) 

ST0R(4,IY) 

ST0R(5,IY) 

ST0R(6,IY) 

ST0R(7,IY) 


TABLE  NUMBER 

NUMBER  OF  SUBLAYER  OF  EACH  LAYER,  NX(IX)  =  THK(IX)/ 
SFRQ(IX) 

ATMOSPHERIC  PRESSURE  IN  MILLIBAR(MB)  USED  IN  SOLVING 
INSOLATION 

BEGINNING  TIME  OF  OUTPUT-TOTAL  NUMBER  OF  HOURS  MINUS  24 
USED  IN  PRINT-OUTPUT  MODULE  2 

EMISSIVITY  BENEATH  AIR  SPACE 

FK(IX)/ALPH(IX)  IN  CAL/CM**2-K 

RICHARDSON  INDEX  NUMBER  USE  IN  SOLVING  CONVECTION 
ENERGY  LOSS 

SURFACE  BENEATH  AIR  SPACE  GEOMETRIC  SHAPE  IN  FRACTION 

(0. 0-1.0) 

RR ( IX) =DELT/ SFRQ**2 .  (FART  OF  HEAT  FLOW  EQUATION) 

SOLAR  AZIMUTH  IN  RADIANS,  SAZ-ATAN  (-COS (DECL)*SIN (TIMER)/ 
( COS (ELF*SIN (DECL) -S IN ( ELF) *COS (TIMER) ) ) ) 

VERTICAL  GRID  SPACING  IN  CM  IN  EACH  LAYER  IX  IN  CM**2/MIN 

INSOLATION  ADJUSTMENT  DUE  TO  ZENITH  ANGLE,  SURFACE  SLOPE 
AND  SURFACE  ASPECT  ANGLE,  SICF-COS (Z) *COS (SLOPE) *SIN (Z) * 

SURFACE  SLOPE  IN  DEGREES  WITH  HORIZONTAL-O  DEGREE,  USED 
IN  SOLVING  INSOLATION 

ABSORPTIVITY  OF  SURFACE  MATERIAL 

EQUALS  SOLAR  INSOLATION  AT  1200  HOURS,  USED  TO  DETERMINE 
IF  SOLAR  INSOLATION  IS  INPUT 

WING  SPEED  IN  CM/SEC 

ESTIMATE  SUB-LAYER  TEMPERATURE  IN  KELVIN 

FK:  HEAT  CONDUCTIVITY  OF  SUB-LAYER  IY  IN  CAL/MIN-CM-K 

RHOC,  FK/ALPH  IN  CAL/CM**2-K 

CONSTANT  DIMENSIONLESS 

INITIAL  SOIL  TEMPERATURE  IN  KELVIN  OF  INITIAL  SOIL  PROFILE 
SAME  AS  STOR(2, IY) 

SAME  AS  STOR(3,IY) 

(Continued) 
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Table  A1  (Continued) 


\ 


V 


Variable 

SUM 

SURFACAZ 

T 

TA 

TAC 

TAK 

TB 

TFREQ 

THK(IX) 

TIME 

TIMER 

TOTTIM 

TPRNT 

TR 

TSK 

TYME 

WATER 

WET 

XL 

XXX(J,1) 

XXX(J,2) 

XXX(J,6) 

XXX(J,3) 

XXX(J,4) 

XXX(J,5) 


Definition 


CALCULATED  INSOLATION  VALUE 

SURFACE  AZIMUTH  IN  DEGREE  WITH  SOUTH  -  0  DEGREE,  USED  IN 
SOLVING  INSOLATION 

SAME  AS  TIME 

AIR  TEMPERATURE  IN  KELVIN 

AIR  TEMPERATURE  IN  DEGREE  CELSIUS 

AIR  TEMPERATURE  IN  KELVIN 

THERMAL  CONDUCTIVITY  OF  BOTTOM  MATERIAL  CAL/CM**2-DEG 
C-MIN 

TIME  STEP  IN  MINUTES  USED  IN  SOLVING  HEAT  FLOW 
LAYER  THICKNESS  IN  CM  OF  LAYER  IX 

TIME  IN  HOURS  IN  WHICH  MATERIAL  TEMPERATURES  ARE  ESTIMATED 
SUN’S  HOUR  ANGLE  IN  RADIANS 

TOTAL  NUMBER  OF  24  HOUR  REPETITIONS  USED  IN  SOLVING 
HEAT  FLOW 

OUTPUT  TIME  PRINT  FREQUENCY  IN  MINUTES 

TEMPERATURE  OF  AIR  SPACE  BENEATH  BOTTOM  MATERIAL 

MATERIAL  SUB-LAYER  TEMPERATURE  IN  KELVIN 

TIME  IN  HOURS  USE  INSOLATION  CALCULATION 

THE  AMOUNT  OF  PRECIPITAL  WATER  IN  MILLIMETRES  (MM) 
CALCULATED  FOR  USE  IN  SOLVING  INSOLATION 

MOISTURE  CONTENT  OF  SURFACE  MATERIAL 

LATENT  HEAT  OF  EVAPORATION  AS  FUNCTION  OF  AIR  AND  GROUND 
TEMPERATURE 

TIME  IN  HOURS  (AIR  TEMPERATURE) 

TIME  IN  HOURS  (RELATIVE  HUMIDITY) 

TIME  IN  HOURS  (WIND  SPEED) 

TIME  IN  HOURS  (AMOUNT  OF  CLOUD  COVER) 

TIME  IN  HOURS 

DEPTH  IN  CENTIMETRES  OF  POINTS  IN  INITIAL  TEMPERATURE 
PROFILE 


(Continued) 
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Wife 


All 


Table  A1  (Concluded) 


Varlabl e 
YYY(J,1) 
YYY(J,2) 

YYY(J,6) 

YYY(J,3) 

YYY(J,4) 

YYY(J,5) 

Z 

ZA 

ZZA 

ZZB 


Definition 


AIR  TEMPERATURE  IN  DEGREES  CELSIUS  TABLE  1 

RELATIVE  HUMIDITY  INPUT  FRACTION,  USED  IN  TABLE  2  SOLVING 
INFRARED  EMISSIONS,  (ATERM)  EVAPORATIVE  HEAT  LOSS  (DTERM) 

WIND  SPEED  INPUT  IN  METRES/ SECOND  AND  CONVERTED  TO 
CENTIMETRE/ SECOND  TABLE  6 

AMOUNT  OF  CLOUD  COVER  IN  FRACTION  (0  TO  1)  USED  IN 
SOLVING  INSOLATION  TABLE  3  INFRARED  EMISSION  (ATERM) 

INSOLATION  IN  CAL/CM**2-MIN,  IF  0.0  AT  12000  HOURS, 
INSOLATION  VALUES  WILL  CALCULATE  TABLE  4 

TEMPERATURE  IN  DEGREES  CELSIUS  AT  POINTS  IN  INITIAL 
TEMPERATURE  PROFILE 

SOLAR  ZENITH  ANGLE,  Z-SIN(DECL)*SIN(ELF)+COS(DECL)* 
COX(ELF)*COS (TIMER) 

SHELTER  HEIGHT  IN  CENTIMETRES  (CM) 

SURFACE  TEMPERATURE  OF  MATERIAL  IN  KELVIN 

BOTTOM  LAYER  TEMPERATURE  OF  MATERIAL  IN  KELVIN 


(Sheet  5  of  5) 


A12 


APPENDIX  B:  TYPICAL  VALUES 
FOR  SYSTEM  DESCRIPTORS 


Table  >1 

toaimcudw  lamut  toittm  for  Seem  Carnot  Hetarlala* 


Subateoca 

Shortwave 

Abeorptivltt 

Longwave 

Sniaeivitr 

Thermal 

Conductivity 

cel/cwJ-mln-*S 

Thtmti 

Dlffutlvlty 

cmJ/nin 

Sandy  eoU 
him 

0.57-0.44 

natural  Hateriala 

0.51-0.55 

0.05 

0.56 

Unfroxam 

— 

— 

0.04 

0.50 

Clayey  eoil 
rroiee 

0.4 

0.45-0.57 

0.15 

0.59 

Unfroaeu 

— 

— 

0.11 

0.24 

Toy  Mil 

Felloe 

0.14-0.95 

O.S 

0.20 

0.510 

Mouthed  field 

0.14-0.9 

— 

— 

— 

Sandy  (ravel 

0.24 

0.24 

0.547 

0.44 

Granite 

0.55-0.45 

0.45 

0.454 

0.96 

Uaeaeoae 

0.7-0.4 

0.5 

0.15 

0.29 

Senditone,  quart i 

0.5-0. 4 

0.54 

0.57-0.72 

0.74 

Freeh  anew 

0.05-0.25 

0.12 

0.05 

0.04 

Old  non 

Clean 

0.55-0.50 

0.54 

0.04 

0.24 

Dirty 

0.50-0.40 

0.54 

— 

— 

let 

0.51 

0.45-0.50 

0.51 

0.70 

Voreet  dull 

O.P2-0.S5 

0.7 

0.010 

0.74 

Concrete  (dry) 

lariated 

ComatructlM  Hateriala 

0.7  0.S4-O.S7 

0.011 

0.17 

Denae 

— 

— 

0.22 

0.43 

Stick 

Sad  Naaonry 

0.55-0.7 

0.5 

0.074-0.12 

0.21-0.35 

Fireclay  (nediuei  brow) 

0.27 

0.75 

0.14 

0.41 

Umber 

Hardwood 

0.4 

0.50 

0.05 

0.07 

Softwood 

— 

— 

0.012-0.022 

0.11 

Tarpaper 

0.55 

0.55 

— 

— 

dabeetoe  abeeta 

O.S 

0.54 

0.017 

0.09 

Aaphalt 

0.4-0.55 

0.54 

0.10 

0.22 

Stalnleaa  eteel 

0.5 

0.12 

5.0 

3.2 

White  planter 

0.07 

0.51 

0.11 

0.2 

Galvaniied  iroa 
lri(bt 

0.45 

0.15 

- 

mm 

Oxidited  (ray 

O.S 

0-24 

— 

— 

Cant  iron 

0.45 

0.44 

4.4 

7.9 

Wrought  iron 

— 

— 

4.4 

9.3 

Glaae 

— 

0.44-0.54 

0.15-0.20 

0.15-0.73 

Folyatyrene 

- 

~ 

0.004 

0.9 

Aluminum 

0.15 

0.04-0.05 

32.0 

3.4 

Solid  HgCOj 

0.04 

Surface  Coetiaae 

0.75 

White  CaO 

0.15 

0.54 

— 

— 

Gray  CaO 

0.57 

0.S7 

— 

— 

White  paiat  (0.04  ca  oa  41) 

0.2 

0.51 

- 

— 

Slack  paint  (0,04  cm  oa  41) 

0.54 

0.44 

- 

- 

Aluminum  paint 

0.2 

0.4 

— 

— 

Clear  varnlah  on  41 

O.S 

0.2 

— 

— 

Suet 

0.54 

0,40-0.54 

— 

— 

*  iiik  1574|  Oka  1974j  Wee haler  and  Uwi  1544j  iMtiM  end  tin  1P4S|  JwOkla  1577i  Aaerlea*  Society  of  kart¬ 
ing,  WiltnitUi,  ut  Air- Condition  Ug  b|UMfi  1577(  uf  peteeael  n— ilcitlw,  0.  Nek,  hum 
national  Foreetry  iMtltuu,  Cniftn  Foreetry  Service,  Chalk  liver,  Ontario,  Canada. 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Thermal  modeling  of  terrain  surface  elements  :  final 
report  /  by  L.  K.  Balick  (Lee  K. )  ...  [et  al.] 
(Environmental  Laboratory,  U.S.  Army  Engineer 
Waterways  Experiment  Station)  i  prepared  for 
Headquarters,  Department  of  the  Any.  —  Vicksburg, 
Miss.  :  U.S.  Army  Engineer  Waterways  Experiment 
Station  ;  Springfield,  Va.  :  available  from  HTI8, 

1981. 

1*2,  [1*3]  p.  :  ill.  1  27  cm.  —  (Technical  report  /  U.8. 
Army  Engineer  Waterways  Experiment  Station  ;  EL-fll-2) 
Cover  title. 

"March  1981." 

"Under  Project  Ho.  L762730ATU2,  Task  Ab,  Work  Unit 
003  and  l»T62719ATb0,  Task  CO,  Work  Unit  006." 
Bibliography:  p.  Ul-1*2. 

1.  Infrared  technology.  2.  Mathematical  models. 

3.  Military  reconnaissance.  It.  Relief  models. 

5.  Temperature.  I.  Balick,  L.  K.  (Lee  K. )  II.  United 
States.  Dept,  of  the  Army.  III.  United  States.  Army 


Thermal  modeling  of  terrain  surface  elements  :  ...  1981. 

(Card  2) 

Engineer  Waterways  Experiment  Station.  Environmental 
Laboratory.  IV.  Title  V.  8eries:  Technical  report 
(United  States.  Army  Engineer  Waterways  Experiment 
Station)  1  EL-Sl-2. 
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